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ABSTRACT 


The  dynamics  of  the  eddy  shedding  by  the  Loop  Current  in  the  Gulf  of  Mexico 
have  been  investigated  using  three  nonlinear  numerical  models:  two-layer, 
barotropic,  and  reduced  gravity.  The  barotropic  and  reduced  gravity  models 
demonstrate  the  individual  behavior  of  the  external  and  internal  modes,  and  provide 
insight  into  how  they  interact  in  the  two-layer  model.  Because  of  the  economy  of 
the  semi- implicit  free  surface  models,  it  was  possible  to  perform  over  100 
experiments  to  investigate  the  stability  properties  of  the  Loop  Current.  Typically 
the  models  were  integrated  3  to  5  years  to  statistical  equilibrium  on  a  1600  x  900 
km  rectangular  domain  with  a  resolution  of  20  x  18.75  km.  Prescribed  inflow 
through  the  model  Yucatan  Channel  was  compensated  by  outflow  through  the  Florida 
Straits. 

A  long-standing  hypothesis  is  that  the  loop  Current  sheds  eddies  in 
response  to  quasi-annual  variations  in  the  inflow.  We  find  that  the  loop 
Current  can  penetrate  into  the  Gulf,  band  westward,  and  shed  realistic 
anticyclonic  eddies  at  almost  an  annual  frequency  with  no  time  variation 
in  the  inflow.  In  this  regime,  the  eddy  shedding  rate  depends  on  the  internal 
Rossby  wave  speed,  an  eddy  diameter  derived  from  conservation  of  potential 
vorticity  on  a  8-plane,  the  angle  of  the  inflow,  and  to  a  lesser  extent  on  the 
Reynolds  number .  Eddy  shedding  can  be  prevented  by  reducing  the  Reynolds  number 
sufficiently.  However,  the  Loop  Current  still  spreads  far  westward.  The  steady 
state  solution  for  a  highly  viscous  case  was  found  to  be  almost  the  same  as  the 
mean  over  an  eddy  cycle  for  a  lower  viscosity  case  which  shed  discrete  eddies  of 

large  amplitude.  Eddy  shedding  and  westward  spreading  of  the  Loop  can  be  prevented 

2 

at  higher  Reynolds  numbers  when  the  beta  Rossby  number,  Rg=vc/ ( BLp  )>2,  where  the 
appropriate  length  scale,  Lp,  is  one-half  the  port  separation  distance  and  vc  is 
the  velocity  at  the  core  of  the  current.  Differential  rotation  (6)  is  also  of 


great  importance  in  determining  the  diameter  and  westward  speed  of  the  eddies  and 
the  penetration  of  the  Loop  Current  into  the  Gulf.  In  a  few  of  the  two-layer 
experiments,  baroclinic  and  mixed  instabilities  were  encountered,  but  experiments 
dominated  by  a  horizontal  shear  instability  of  the  internal  mode  produced  the  most 
realistic  results.  For  sufficiently  high  Reynolds  numbers  the  shear  instability 
occurred  in  both  the  barotropic  and  reduced  gravity  models.  However,  for  realistic 
parameter  values  eddy  shedding  occurred  in  the  two-layer  and  reduced  gravity 
models,  but  not  in  the  barotropic  model. 

Consistent  with  potential  vorticity  conservation  dynamics,  the  Loop  Current 
and  its  eddy  shedding  behavior  was  quite  insensitive  to  the  location  and  width  of 
the  inflow  and  outflow  ports,  so  long  as  the  western  boundary  did  not  interfere 
with  the  shedding  process  and  the  ports  were  not  separated  by  much  less  than  1/V2” 
times  the  theoretical  eddy  diameter,  i.e.,  when  RB  <  2.  When  the  entire  eastern 
boundary  was  opened,  the  outflow  remained  confined  to  a  current  adjacent  to  the 
southern  boundary.  Also,  while  the  solution  depends  on  the  maximum  velocity  at 
inflow,  it  is  relatively  insensitive  to  the  shape  of  the  inflow  profile. 

In  the  presence  of  significant  deep-water  inflow  through  the  Yucatan  Straits, 
bottom  topography  may  prevent  Loop  Current  penetration,  westward  spreading,  and 
eddy  shedding.  In  these  cases  the  interaction  between  the  bottom  topography  and 
the  pressure  field  near  the  Florida  Shelf  results  in  a  near  balance  between  the 
pressure  torques  and  the  nonlinear  terms  in  the  mass  transport  vorticity  equation. 
When  the  Yucatan  Straits  deep  water  inflow  is  reduced  or  the  Florida  Shelf  is  moved 
to  the  east,  the  eddy  shedding  reappears.  A  kinematic  analysis  shows  that  a 
sufficiently  strong  current  following  f/h  contours  of  the  Florida  Shelf  and 
intersecting  the  Loop  Current  at  large  angles  can  locally  prevent  northward 
penetration  of  the  Loop  Current  and  effectively  reduce  the  port  separation.  Thus, 
the  effect  of  the  Florida  Shelf  is  similar  to  cases  in  the  reduced  gravity  model 
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where  the  ports  are  too  close  for  eddy  shedding  to  occur,  i.e. ,  when  %  >2. 
Bottom  topography  also  inhibits  development  of  baroclinic  instability,  yielding 
solutions  more  closely  resembling  those  from  the  reduced  gravity  model  than  from 
the  two-layer  flat  bottom  model.  However,  movement  of  the  shed  eddies  is 
significantly  modified  by  the  introduction  of  topography. 

In  the  presence  of  realistic  time  variations  in  the  upper  layer  inflow,  the 
eddy  shedding  period  is  dominated  by  the  natural  period  rather  than  the  forcing 
period,  although  the  influence  of  the  latter  is  not  negligible. 
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1.  Introduction 


The  Gulf  of  Mexico  is  remarkable  for  the  variability  and  intensity  of  its 
circulation.  The  Loop  Current  in  the  eastern  Gulf  is  of  particular  interest  since 
it  supplies  a  significant  fraction  of  the  source  water  for  the  Gulf  Stream  - 
approximately  30  x  10^  m3 sec- 1  (Nowlin,  1972) .  Entering  the  Gulf  through  the 
Yucatan  Straits,  the  Loop  Current  traces  an  anticyclonic  path  which  may  extend  far 
to  the  north  before  turning  southward  and  exiting  through  the  Straits  of  Florida. 
Large  anticyclonic  eddies  have  been  observed  to  break  off  from  the  Loop  Current 
following  these  intrusions  (Fig.  1) .  The  eddies  typically  have  diameters  of  about 
360  km  and  translate  into  the  western  Gulf  with  a  mean  speed  of  2.4  an  sec  ■*" 
(Elliott,  1979) .  One  hypothesis  holds  that  the  persistent  anticyclonic  gyre  in  the 
western  Gulf  is  periodically  renewed  by  these  break-off  eddies  (Cochrane,  1972; 
Behringer,  et  al.  1977).  This  contrasts  with  the  suggestion  by  Sturges  and  Blaha 
(1976)  that  the  western  gyre  may  be  driven  by  wind-stress  curl. 


Cochrane  (1965)  proposed  that  the  path  for  the  loop  Current  might  vary  in 
response  to  an  annual  variation  in  the  strength  and  cross-stream  structure  of  the 
surface  currents  in  the  Yucatan  Straits.  Leipper  (1970) ,  Cochrane  (1972) ,  Maul 
(1977) ,  and  others  have  contributed  limited  observational  evidence  for  a  seasonal 
cycle  in  the  Loop  Current,  although  each  author  identified  considerable 
year-to-year  variability  in  the  data.  Recently,  Molinari  (1978)  has  challenged  the 
hypothesis  of  an  annual  cycle  in  the  Loop  Current: 


Circulation  patterns  in  the  Gulf  given  by  climatological 
averages  and  several  synoptic  studies  suggest  an  annual 
cycle  for  the  intrusion  of  the  Loop  Current.  The  Loop 
intrudes  north,  in  these  representations,  into  the  Gulf 
during  late  winter  and  spring,  reaches  maximum  pene¬ 
tration  in  the  Gulf  during  the  summer,  sheds  an  eddy 
which  drifts  west  during  the  late  sumner  and  fall  and 
exhibits  minimum  penetration  in  the  winter.  A  reanalysis 
of  the  historical  data  set  suggests  that  the  climatological 
results  are  severely  biased  by  the  temporal  sampling 
distribution.  In  fact,  large  sunmer-fall  intrusions 


occurred  in  1966 ,  1969,  1973,  and  1974,  followed  by  maxinun 
penetration  in  the  winter,  rather  than  simmer. 


Glaring  inadequacies  both  in  the  data-base  outside  the  Loop 
Current  and  our  knowledge  of  the  dynamics  of  the  circulation 
in  the  Gulf  are  apparent.  Circulation  patterns  in  the  absence 
of  the  Loop  or  detached  eddies  are  unknown.  Thus,  the  magnitude 
and  region  of  exchange  between  the  eastern  and  western  Gulf  are 
unknown.  Little  data  (even  less  current  meter  data)  exist  to 
define  the  deep  circulation.  The  effect  of  cold  fronts  on  the 
circulation,  the  effect  of  the  boundary  conditions  at  the  Yucatan 
and  Florida  Straits  on  the  Loop  Current,  and  the  effect  of  the 
wind  on  the  circulation  have  yet  to  be  quantified.  As  basic  a 
question  as  what  drives  the  western  Gulf  circulation,  detached 
eddies  or  curl  of  the  wind  stress  is  still  unanswered.  The  time 
appears  appropriate  for  the  use  of  simple  analytical  and/or 
numerical  models  to  address  these  questions  and  to  provide  the 
basis  for  more  process  oriented  observational  studies  rather 
them  the  circulation  mapping  exercises  of  the  past. 

The  purpose  of  this  paper  is  to  simulate  and  understand  the  dynamics  of  the 
Loop  Current  and  its  eddy  shedding  using  an  ensemble  of  numerical  models,  including 
two-layer,  barotropic,  and  reduced  gravity  prototypes.  A  series  of  inportant 
inter-related  questions  are  addressed:  Can  a  simple  model  simulate  the  Loop  Current 
and  eddy  shedding  in  statistical  agreement  with  the  limited  observations?  Is  the 
eddy  shedding  rate  controlled  by  internal  dynamics  or  by  external  changes  in 
boundary  conditions?  What  determines  the  scale  of  the  shed  eddies?  What  dynamics 
control  the  Loop  Current  penetration  and  eddy  break-off?  How  does  the  magnitude 
and  cross-stream  profile  of  the  current  through  the  Yucatan  Straits  affect  Loop 
Current  behavior?  What  happens  to  the  eddies  once  they  are  shed?  What  is  the 
influence  of  bottom  topography  on  eddy  shedding  and  the  Loop  Current? 


Because  of  the  efficiency  of  the  primitive-equation  model  formulation,  in 
which  free-surface  and  internal  gravity  waves  are  treated  implicitly,  over  100 

model  experiments  covering  a  wide  range  of  the  parameter  space  ware  integrated 
to  statistical  equilibrium.  Previous  studies  using  a  barotropic  model  (Paskausky 

and  Reid,  1972)  and  a  quasi-geostrophic  two-layer  model  (Wert  and  Reid,  1972) 
were  extremely  limited  in  the  range  of  parameter  space  studied  and  the  models  were 
not  integrated  to  statistical  equilibrium. 
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The  formulation  and  parameters  of  the  three  models  used  in  this  study  are 
presented  in  section  2.  Numerical  aspects  of  model  design  are  discussed  in 
Appendix  B,  where  particular  attention  is  given  to  the  semi-implicit  treatment  of 
the  gravity  waves,  to  the  boundary  conditions  for  the  resulting  Helmholtz 
equations,  and  to  the  boundary  conditions  at  the  open  ports.  In  section  3  we  first 
explore  the  capability  of  the  nixnerical  models  to  perform  realistic  simulations  of 
the  Loop  Current  and  the  eddy  shedding.  Using  realistic  parameters,  we  determine, 
which,  if  any,  of  the  models  produce  eddies  with  realistic  size,  amplitude,  and 
movement  and  the  observed  quasi-annual  period  of  eddy  shedding. 

In  section  3b  we  vary  the  model  parameters  to  assess  the  influence  of  the 
internal  dynamics  on  the  eddy  shedding  by  the  Loop  Current  and  on  the  character  of 
any  shed  eddies.  Regression  analysis  is  used  to  test  a  number  of  dynamical 
hypotheses  based  on  Rossby  wave  theory  and  on  conservation  of  potential  vorticity 
on  a  0-plane .  In  section  3c  we  compare  a  steady  solution  with  a  low  Reynolds 
number  to  a  linear  viscous  solution  and  to  the  mean  of  a  case  which  exhibited 
eddy  shedding.  The  regimes  found  in  a  barotropic  flat  bottom  model  are  discussed 
in  section  3d  and  compared  to  those  of  the  reduced  gravity  model,  noting  the  role 
of  the  divergence  term  in  the  vorticity  equation  in  creating  the  differences.  In 
section  3e  we  compare  the  relative  realism  of  results  from  cases  with  horizontal 
shear  instability,  baroclinic  instability,  and  mixed  stability. 

In  section  4  we  investigate  the  importance  of  external  influences  on  the  eddy 
shedding  by  the  Loop  Current.  In  the  three  subsections  we  explore  the  response  of 
the  model  (a)  to  changes  in  the  location,  width,  and  separation  of  the  inflow  and 
outflow  ports,  (b)  to  the  introduction  of  bottom  topography,  and  (c)  to  spatial  and 
temporal  changes  in  the  inflow  transport.  In  section  4a  some  of  the  results  for 
the  reduced  gravity  and  barotropic  flat  bottom  models  are  summarized  in  a  stability 
regime  diagram.  In  section  4b  a  kinematic  analysis  is  used  in  explaining  the 
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important  role  of  the  Florida  Shelf  when  there  is  sufficient  deep  water  inflow 
through  the  Yucatan  Straits.  The  influence  of  the  topography  on  the  baroclinic 
stability  of  the  Loop  Current  is  also  considered. 


2.  The  Gulf  of  Mexico  models 


Three  models  are  used  in  this  study:  two-layer,  barotropic,  and  reduced 
gravity.  The  last  two  are  mathematically  identical  but  have  differing 
gravitational  accelerations.  In  the  reduced  gravity  model  a  value  appropriate  for 
the  first  internal  mode  is  used  and  the  lower  layer  is  assumed  infinitely  deep  and 
at  rest.  The  barotropic  and  reduced  gravity  models  demonstrate  the  individual 
behavior  of  the  external  and  internal  modes  and  provide  insight  into  how  they 
interact  in  the  two-layer  model.  They  also  allow  us  to  investigate  certain 
phenomena  in  the  simplest  context. 

The  two-layer  and  reduced  gravity  models  are  assumed  to  be  stably  stratified 
and  a  fixed  density  contrast  is  specified  between  the  two  immiscible  layers.  All 
other  thermodynamic  effects  are  neglected.  For  all  three  models  we  assume  a 
hydrostatic,  Boussinesq  fluid  with  a  free  surface  in  a  rotating  right-handed 
coordinate  system.  Retaining  only  the  first  term  in  the  perturbation  expansion  for 
the  advective  terms,  the  vertically  integrated  equations  of  motion  used  in  the 


models  are 


r  +  <’-?i +  -v>7f * f  ?i  =  -V  pi 


+  Ti  '  Ti+1  +  A  V2 


9hi  +  7.^-0  , 
3 1  1 


where  i  *  1,  2  for  the  two- layer  model,  i=l  for  the  barotropic  and  reduced  gravity 


models  and 


v  »=  jL  i  +  JL  ^ 

3x  3y  J  » 


P]  =  9h1  ; 

P2  ■  P!  -  g‘hT  » 

✓v  A 

^  a  h1vi  =  hi  (u,  1  +  vi  j). 
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9  ^2  '  pl) 

P  * 

f  •  fo  ♦  e(y  -  J0)  , 

Ti  -  Ti  ’  ♦  Tiy  j 

See  Appendix  A  for  symbol  definitions  . 

In  the  reduced  gravity  model,  the  lower  layer  momentum  equation  degenerates  to 
VP2=0,  and  in  the  barotropic  model  t2  becomes  the  bottom  stress.  Laplacian 
horizontal  friction  is  such  a  crude  representation  of  sub-grid  scale  motion  that  we 
did  not  believe  it  necessary  to  use  the  more  complicated  form  which  actually 
results  from  the  derivation  of  the  model  equations.  However,  a  number  of 
experiments  were  performed  in  which  AH72v  rather  than  AV2V  was  used  for  horizontal 
friction.  Only  small  quantitative  differences  were  found  except  in  the  two-layer 
model  with  bottom  topography  and  significant  deep  water  transport.  There  the 
regime  boundary  between  steady  and  eddy  shedding  regimes  was  shifted.  For  further 
details,  see  Blake,  et  al.  (1980) . 

Figure  2  shows  the  model  data  in  superimposed  on  a  map  of  the  Gulf  of  Mexico . 

We  use  a  1600  x  900  km  rectangular  domain,  neglecting  the  20°  counter-clockwise 
rotation.  Thus,  the  Coriolis  parameter,  f=f (y)  only.  We  use  the  &-plane 
approximation  for  that  variation.  The  model  is  driven  from  rest  by  inflow  through 
the  southern  port  (Yucatan  Straits) .  This  is  exactly  compensated  by  outflow 
through  the  eastern  port  (Florida  Straits) .  In  most  experiments  the  inflow  port  is 
160  km  wide,  centered  400  km  from  the  eastern  boundary.  The  outflow  port  is  150  km 
wide,  centered  75  km  from  the  southern  boundary.  In  a  few  cases  the  Florida  and 
Yucatan  Straits  are  explicitly  included  as  extensions  to  the  model  domain.  Wind 
forcing  is  neglected  in  the  experiments  described  here.  This  does  not  imply  that 
we  belive  wind  forcing  is  uninportant  to  the  Gulf  circulation.  Rather,  we  have 
chosen  the  Loop  Current  forcing  as  the  focus  of  this  study.  Figure  3  shows  the 
idealized  bottom  topography  of  the  Gulf  used  in  some  experiments. 

_ 9 _  . 


Except  at  the  inflow  and  outflow  ports,  the  walls  are  rigid  and  the  no-slip 


condition  is  prescribed  on  the  tangential  flow.  At  inflow  the  profile  of 
transport,  V^,  is  prescribed.  In  most  cases  this  profile  is  parabolic.  Because 
the  current  is  nearly  geostrophic,  the  upper  layer  velocity  maximum  is  west  of  the 
center  of  the  Yucatan  Straits.  The  normal  flow  at  the  eastern  port  is 
self-aetermined  with  the  integral  constraint  that  the  total  transport  out  from  each 
layer  match  the  inflow.  Usually  the  amplitude  of  the  tangential  component  was  set 
to  zero  one-half  grid  distance  outside  the  physical  domain  at  both  ports.  This 
weak  over-specification  at  outflow  eliminates  the  possibility  that  the  Loop  Current 
will  exit  the  basin  at  unrealistic  angles.  More  realistically,  outflow  should 
occur  through  a  channel  modeling  the  Florida  Straits.  This  was  done  in  a  few 
cases.  See  Appendix  B  for  further  discussion  of  the  design  of  the  semi- implicit 
numerical  models. 

The  parameters  of  a  'Standard"  case  for  each  model  are  given  in  Table  1.  The 
value  of  g'  is  the  equivalent  for  a  reduced  gravity  model.  In  the  two- layer  model 
this  value  is  multiplied  by  (H^+H2)/H2  to  yield  the  same  internal  values  for  the  gravity 
wave  speed,  radius  of  deformation,  and  Rossby  wave  speed  as  in  a  reduced  gravity 
model.  The  inflow  transport  is  spun  up  with  a  time  constant  of  30  days.  The  only 
dissipation  in  the  model  is  due  to  horizontal  friction.  Interfacial  and  bottom 
stresses  were  excluded  in  almost  all  cases.  Two  cases  were  selected  in  which  the 
potential  role  of  bottom  friction  was  the  greatest.  These  were  repeated  including 
quadratic  bottom  friction  and  the  results  are  discussed  in  section  3.  Our  standard 
case  value  for  the  horizontal  eddy  viscosity  is  larger  than  required  for  stable 
integration  of  the  model.  This  value  yielded  a  constant  eddy  shedding  period  for 
the  Loop  Current.  Smaller  values  produced  some  irregularity  in  this  period  without 
substantially  altering  the  long  term  statistical  properties  of  the  eddy  shedding 
process.  Although  lower  eddy  viscosities  were  used  in  some  cases,  the  larger  value 
was  used  in  most  to  decrease  the  length  of  the  integration  required  to  obtain 
stable  statistics  and  to  facilitate  the  analysis  of  the  results. 
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Table 

1.  Model  Parameters 

for  Standard 

Case 

A 

lO^cm^sec-! 

6 

2  xl0-£ 3cm~£ sec“£ 

£0 

5  x  10~^sec_£ 

P 

1  gm  cm-^ 

9 

980  cm  sec"2 

-* 

Ti 

0 

9' 

3  cm  sec“2 

Ax 

20  km* 

H1 

200  m 

Ay 

18.75  km* 

h2 

2800  m 

At 

1 . 5  hr 

Domain 

Size,  xL  by  yL 

1600  x 

900  km 

Southern  Port  Width,  L  160  km 

Eastern  Port  Width,  Le  150  km 

Center  of  southern  port  at  Xp  1200  km 

Center  of  eastern  port  at  yp  75  km 

Upper  Layer  Inflow  Transport**  20  x  lO^m^sec-^  (20  Sv) 

Lower  Layer  Inflow  Transport  10  x  lO^m^sec-^  (10  Sv) 

Angle  of  inflow  from  x-axis,  9^  90° 

For  the  barotropic  model  the  initial  maximum 
depth  is  H=3000  m  and  the  inflow  transport  is  30  Sv. 

*  for  a  given  variable 

**also  for  the  standard  reduced  gravity  model 


3.  The  influence  of  internal  dynamics  on  eddy  shedding  by  the  Loop  Current 

Our  first  goal  was  to  simulate  the  Loop  Current  and  its  eddy  shedding  in  a 
manner  which  resembles  the  observed  behavior.  A  nonlinear  two-layer  hydrodynamic 
model  was  chosen  as  the  sinplest  model  likely  to  succeed.  A  second  active  layer 
allowed  the  possibility  of  baroclinic  instability  and  other  modal  interactions,  and 
allowed  inclusion  of  both  an  idealized  pycnocline  and  idealized  topography  of  the  Gulf. 

Once  we  demonstrated  the  ability  of  the  model  to  perform  a  realistic 
simulation,  we  1)  examined  the  dynamics  of  the  eddy  shedding  Loop  Current  in  the 
simplest  context,  2)  looked  for  the  existence  of  other  regimes  within  the 
neighboring  parameter  space,  and  3)  examined  the  effects  of  topography,  time 
variations  in  the  inflow,  and  changes  in  the  port  configuration.  The  third  area  of 
investigation  is  discussed  in  section  4. 

a.  A  simulation  of  the  eddy  shedding  by  the  Loop  Current 

A  longstanding  hypothesis  is  that  the  Loop  Current  sheds  eddies  in  response  to 
quasi-annual  variations  in  the  vorticity  and  the  inertial  character  of  the  current 
through  the  Yucatan  Straits.  These  variations  in  the  total  transport,  or  the 
horizontal  and  vertical  structure  presumably  affect  the  penetration  distance  of  the 
Loop  Current.  When  the  current  retreats,  an  eddy  is  left  behind.  Despite  the 
popularity  of  this  idea,  we  first  searched  for  purely  internal  mechanisms  where  the 
eddy  shedding  results  from  unstable  configurations  of  the  Loop  Current,  independent 
of  any  finite  amplitude  perturbation  due  to  time  variations  in  the  inflow. 

Instead  of  starting  with  the  sinqplest  possible  model  and  progressing  toward 
more  complicated  experiments,  we  designed  our  first  experiment  to  be  as  realistic 
as  possible  within  the  context  of  the  two-layer  model  and  a  time  invariant  inflow 
profile  for  V.  We  then  looked  for  the  simplest  model  which  produced  essentially 
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the  same  results.  Hie  parameters  for  the  first  ramerical  experiment  using  the  two- 
layer  model  are  listed  in  Table  1,  except  that  the  upper  layer  inflow  transport  is  25 
Sv,  and  the  lower 'layer  5  Sv.  The  idealized  topography  for  the  Gulf  of  Mexico  shown 
in  Figure  3  is  also  included.  The  inflow  transport  is  steady  except  for  spin-up 
with  a  time  constant  of  30  days. 

Figure  4  shows  a  sequence  of  snapshots  of  the  pycnocline  anomaly  at  70-day 
intervals.  The  pycnocline  anomaly  is  the  deviation  of  the  interface  from  the 
initial  horizontal  surface.  Downward  deviations  are  positive  (upper  layer 
thickness  greater  than  initial) .  Figure  4a  shows  the  Loop  Current  already 
penetrating  far  into  the  Gulf  and  beginning  to  develop  an  eddy  structure.  Figure 

4b  shows  an  eddy  starting  to  break  off  from  the  loop.  Figure  4c  shows  the  Loop 

Current  and  eddy  just  after  separation.  At  this  stage  the  Loop  Current  is  a 

source-sink  flow  with  less  penetration  into  the  Gulf.  Figure  4d  shows  the  eddy  has 

moved  westward  and  the  Loop  Current  has  penetrated  further  into  the  Gulf. 

Returning  to  Figures  4a  and  4b,  we  see  that  when  an  eddy  reaches  the  western 
boundary,  it  moves  northward  with  final  decay  in  the  northwestern  part  of  the 
basin.  Figure  5a  shows  the  pycnocline  as  a  function  of  time  at  the  longitude 
x  =  1010  km  from  the  western  boundary. 

Even  though  the  inflow  profile  for  V  is  steady,  the  eddy  shedding  occurs  in  a 
regular  manner  with  a  period  of  290  days.  This  period  compares  favorably  with  the 
quasi-annual  period  observed.  Furthermore,  when  parameters  suitable  for  the  Gulf 
of  Mexico  are  used,  the  size,  amplitude,  and  movement  of  the  eddies  are  also 
consistent  with  observations  (Nowlin,  1972;  Cochrane,  1972;  Elliott,  1979).  Of 
particular  interest  is  the  realistic  value  for  the  eddy  shedding  period  predicted 
by  the  model  with  no  time  variation  in  the  inflow,  since  the  popular  hypothesis  is 
that  the  Loop  Current  exhibits  a  quasi-annual  shedding  period  in  response  to  the 
seasonal  variations  in  the  inflow  through  the  Yucatan  Straits. 


Thus  far,  we  have  shown  that  the  model  can  realistically  simulate  observed 
features  of  the  eddy  shedding  by  the  Loop  Current  when  model  parameters  realistic 
for  the  Gulf  of  Mexico  are  used.  In  addition,  we  have  shown  that  time  variations 


in  the  inflow  are  not  required  for  the  eddy  shedding  to  occur  at  a  realistic  rate. 


At  this  point  we  want  to  examine  the  dynamics  of  the  eddy  shedding  by  the  Loop 
Current  in  the  simplest  context.  Figure  5  shows  the  pycnocline  anomaly  vs.  time  at 
a  longitude  190  km  west  of  the  center  of  the  inflow  port  for  three  different  cases. 
Figure  5a  is  for  the  case  with  bottom  topography  already  described.  Figure  5b  is 
for  a  two-layer  flat  bottom  case  using  the  parameters  given  in  Table  1.  Figure  5c 
is  for  a  reduced  gravity  experiment  using  the  same  parameters,  except  that  the 
lower  layer  is  infinitely  deep  and  at  rest.  Most  striking  is  that  all  three 
experiments  produce  nearly  the  same  results  in  terms  of  eddy  shedding  period,  eddy 
diameter,  and  eddy  amplitude,  within  20%  or  less  in  almost  all  respects.  In  all 
three,  discrete  eddies  are  shed  in  a  regular  manner  with  little  evidence  of  other 
oscillations.  The  most  notable  difference  is  that  the  case  with  the  bottom 
topography  of  Figure  3  did  not  begin  to  shed  eddies  for  almost  three  years.  The 
barotropic  model  evolved  to  a  steady  state  without  shedding  eddies  in  this 


parameter  range  with  30  Sv  inflow  and  3000  m  depth  (see  Section  4,  Fig.  23a). 

The  ability  of  the  reduced  gravity  model  to  produce  essentially  the  same 
results  as  the  two-layer  model  is  sufficient  to  demonstrate  that  the  discrete  eddy 

shedding  regime  found  in  these  three  cases  is  not  due  to  baroclinic  instability, 
but  instead  is  due  to  a  horizontal  shear  instability  of  the  internal  mode.  The 
reduced  gravity,  two-layer  flat  bottom,  and  two-layer  model  with  idealized  Gulf  of 
Mexico  topography  do  not  agree  in  all  the  dynamical  regimes  we  found,  as  we  shall 
see  later.  However,  they  do  agree  for  a  regime  in  accord  with  observed  features  of 
the  Loop  Current-eddy  system.  This  suggests  that  the  simplest  useful  model  to 
study  the  dynamics  of  the  eddy  shedding  from  the  Loop  Current  is  the  reduced 
gravity  model. 

Before  turning  to  the  reduced  gravity  model,  we  will  demonstrate  the 
importance  of  differential  rotation  (8)  to  the  eddy  shedding.  Figure  6 
shows  the  free  surface  anomaly  at  day  1080  for  an  f-plane  case  using  the  two-layer 
flat  bottom  model  with  the  parameters  of  Table  1,  except  that  8=0  and  the  inflow 
angle  is  27°  west  of  normal.  Even  with  the  westward  inflow  component,  the  solution 
evolved  to  the  steady  state  source-sink  flow  shown  in  Figure  6  without  shedding  an 
eddy.  The  experiment  was  continued  after  adding  a  time  varying  inf  lew  with  an  annual 
signal  much  larger  than  observed  (total  inflow  transport  ranging  from  20  to  40  Sv) . 
Still  no  eddy  shedding  occurred.  With  a  sufficiently  high  Reynolds  number,  no 
doubt  seme  physical  instability  would  occur.  Hcwever,  it  is  clear  that  8  plays  an 
essential  role  in  the  discrete  eddy  shedding  regime  shown  in  Figures  4  and  5. 

A  large  number  of  numerical  experiments  were  conducted  using  the  reduced 
gravity  model.  All  of  those  with  a  steady  inflow,  time  invariant  values  of  other 
specified  parameters,  and  port  geometry  not  grossly  different  from  the  Gulf  of 
Mexico  are  listed  in  Table  2.  Barotropic  flat-bottom  cases  are  also  included. 
Fortunately,  we  were  able  to  perform  enough  reduced  gravity  experiments  to  test  a 
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Case  Differences  from  Table  1  dm  ye  bm  pm  cm  v£n  vmxe  he  hm  ha  Re  R|)  Notes 


Symbol  Definitions  for  Table  2 


N,n  Stable  regime  for  the  reduced  gravity  model  (N)  and  the  barotropic  model  (n) 
in  which  there  was  no  westward  spreading  of  the  Loop  Current  and  no  eddy 
shedding. 

W,w  Stable  regime  for  the  reduced  gravity  model  (W)  and  the  barotropic  model  (w) 
in  which  there  was  westward  spreading  of  the  Loop  Current,  but  no  eddy 
shedding. 

<^j,  North-south  diameter  in  km  of  an  eddy  measured  between  speed  maxima  as  the 
center  passed  110  km  west  of  the  western  boundary  of  the  inflow  port.  For 
regime  W,w  all  quantitites  were  calculated  treating  the  loop  as  we  would 
an  eddy. 

ye  Distance  from  the  southern  boundary  in  km  of  the  eddy  pressure  centers  as 
they  passed  a  longitude  10  km  west  of  the  western  boundary  of  the  inflow 
port.  In  regime  N,n  quantities  were  calculated  at  whatever  longitude  the 
extreme  value  occurred. 

fcta  Distance  from  the  southern  boundary  in  km  of  the  northern  speed  maximum  as 
the  center  of  an  eddy  passed  110  km  west  of  the  western  boundary  of  the 
inflow  port. 

Eddy  shedding  period  in  days.  In  the  stable  cases  the  type  of  stable  regime 
is  listed. 

qn  Average  westward  speed  in  an/sec  of  an  eddy  between  110  km  and  310  km  west 
of  the  western  boundary  of  the  inflow  port.  In  cases  marked  by  an  asterisk 
(*)  this  was  not  clearly  defined  due  to  interaction  with  other  eddies. 

v.  Maximum  current  speed  at  inflow  in  cm/ sec. 
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Vnixe  Maximum  current  speed  in  cm/ sec  110  km  west  of  the  western  boundary  of  the 
inflow  port  as  an  eddy  passed  through  that  longitude. 

he  Representative  depth  in  m  of  the  interface  in  an  eddy  as  the  center  passed 
110  km  west  of  the  western  boundary  of  the  inflow  port.  This  was  taken  as 
the  mean  of  the  values  at  the  northern  and  southern  speed  maxima. 

hm  Maximum  depth  in  m  of  the  interface  in  an  eddy  as  the  center  passed  a 

longitude  110  km  west  of  the  western  boundary  of  the  inflow  port.  In  the 
barotropic  cases  it  is  the  maximum  FSA  in  cm. 

ha  Amplitude  in  m  of  the  eddy  pycnocline  deformation  as  the  eddy  center  passed 
10  km  west  of  the  western  boundary  of  the  inflow  port.  This  was  measured 
relative  to  the  ambient  value  of  the  pycnocline  anomaly  surrounding  an  eddy 
and  is  not  necessarily  the  same  as  the  maximum  pycnocline  anomaly  in  an 
eddy.  In  the  barotropic  cases  we  used  FSA  in  cm. 

Re  Reynolds  number,  v^h/A  where  Lh  is  half  the  width  of  the  southern  (inflow) 
port. 

2 

Pg  Beta  Rossby  number,  v^/(8L£)  where  Lp  is  half  the  distance  between  the 

centers  of  the  ports.  When  the  entire  eastern  boundary  is  open,  Lp=°°.  In 
cases  marked  by  D2  an  eastern  port  centered  75  km  south  of  the  western 
boundary  of  the  inserted  land  mass  was  used  (Fig.  17) . 

The  values  in  this  table  were  measured  to  within  3%  or  better,  except  for  ha  which 

was  only  determined  within  about  10%. 
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Additional  information  on  differences  from  Table  1 

Units:  A  in  1 0®  am^sec  \  8  in  10  ^  an  ^sec  f  in  10  '‘sec  \  inflow  in  Sv,  L, 

L  ,  x, ,  y_ ,  x  ,  y  in  km 
e'  l  ■'L  p'  Jp 

Lq  =  .9375L  unless  explicitly  given 

Dl  Hie  Yucatan  and  Florida  Straits  were  added  to  the  model  domain  (Fig.  20) . 

In  these  cases  the  value  of  8^  measured  from  the  x-axis  is  an  average  over 
an  eddy  cycle  at  the  entrance  to  the  Gulf;  vin  is  also  measured  at  the 
entrance  to  the  Gulf.  At  the  entrance  to  the  Yucatan  Straits,  the  inflow 
was  in  the  y-direction. 

D2  A  land  area  was  inserted  to  match  the  location  of  the  Florida  Shelf 
(Fig.  17). 


Notes  on  results  in  Table  2 

Irregular  eddy  shedding  period:  322-438  days,  std.  dev.  -  37  days;  std  dev 
for  djjj  *  5  km,  for  ■  13  tan.  Values  in  the  table  are  means  for  the  last 
12  eddy  cycles  in  a  15-yr  integration.  The  range  of  eddy  periods  is  similar 
to  that  observed  (Behringer,  et  al.,  1977). 

Irregular  edd^y  shedding  period:  328-352  days,  std  dev  -  8  days.  The  value 
in  the  table  is  the  mean  over  the  last  5  eddy  cycles  in  a  7-year 
integration. 


These  experiments  did  not  reach  statistical  equilibrium  within  the  5  to  10 
year  integrations  (see  section  4a) .  Values  in  the  table  are  for  the  5th 
eddy  or  eddy  cycle.  Significant  rates  of  change  are  given  in  percent  per 
eddy  cycle  for  each  case: 


HG44: 

RG45: 


BG46: 


-.8  for  <%,  -.3  for  bin,  2.5  for  p^, 

-3  for  d^,  -2  for  bm,  4.5  for  Rn'  1*5  for  Cm, 
-3  for  Vmxe*  -1.5  for  ye, 

-.6  for  dm,  2.7  for  r„,  -1  for  c^,  -.4  for  ye 


Some  values  may  not  be  representative  due  to  the  distortion  of  the  eddies  as 
they  interacted  with  the  strong  secondary  circulation. 
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number  of  dynamical  hypotheses  using  regression  analysis.  In  all  cases  sufficient 
data  is  included  in  the  tables  to  repeat  the  analysis.  Of  the  cases  listed  in 
Table  2,  35  were  used  in  the  regression  analyses.  Only  the  Yucatan  Peninsula  case, 
the  cases  with  the  entire  eastern  boundary  open,  the  stable  cases,  and  the 
baro tropic  cases  were  excluded,  unless  otherwise  noted. 

Once  an  eddy  is  shed,  it  is  observed  to  move  westward.  Using  data  from  the 

tables,  we  found  that  the  linear  correlation  between  the  westward  component  of  the 
eddy  movement  and  the  non-dispersive  internal  Ftossby  wave  speed,  c.^r=Bg'h^/f  ,  is 

0.991.  However,  averaged  over  the  35  cases  used,  this  theoretical  speed  is  40%  greater 
than  observed  in  the  numerical  model,  when  he  from  Table  2  is  used  for  h]_. 

Including  the  dispersive  contribution  to  the  phase  speed  reduced  this  difference  to 
<7%,  but  degraded  the  correlation,  although  it  was  still  -0.90.  Using  he  fran 
Table  2  yielded  better  agreement  than  using  the  mean  thickness  of  the  upper  layer 
or  the  maximum  thickness  at  the  center  of  the  eddy,  indicating  some  nonlinear 
contributioi  to  the  westward  speed  of  the  isolated  eddies.  McWilliams  and  Flierl 
(1979)  discuss  the  propagation  and  dynamics  of  isolated  eddies  in  a  variety  of 
situations,  most  more  nonlinear  than  evolved  in  this  study. 

Reid  (1972)  presents  an  argument  based  on  conservation  of  potential  vorticity 
on  a  0-plane  to  estimate  the  penetration  of  the  Loop  Current  into  the  Gulf  of 
Mexico.  We  have  found  an  extension  of  that  argument  and  its  application  to  be 
useful  in  understanding  a  number  of  dynamical  aspects  of  the  numerical  solutions. 

In  a  frictionless,  steady  state,  reduced  gravity  model 

^1  *  ^  *  constant  (4) 

hl 

along  a  streamline,  where  ^  is  the  relative  vorticity  in  the  active  layer.  If  we 
assume  the  flow  is  geostrophically  balanced,  then  contours  of  hj  are  streamlines, 
and  with  the  0-plane  approximation  (4)  becomes 
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5i +  6  y =  s 


(5) 


where  we  have  set  yQ  =  o  and  c0  is  the  relative  vorticity  when  y  =  o.  Writing  5  in 
natural  coordinates  and  neglecting  variations  in  radius  of  curvature 


C 


1 


(6) 


at  the  core  of  the  current  where  there  is  no  contribution  from  shear,  where  vc  is 
the  speed  at  the  core  of  the  current,  ds  is  an  element  of  arc  length  positive  in  a 
clockwise  direction,  and  9  is  the  angle  between  the  velocity  vector  and  the 
positive  x-axis.  Noting  &  =  S1-n  e>  (5)  becomes 


sin  @  ‘37 


=  '^y 


_co  .  de 

V7s'n  6o  d? 


y=0 


(7) 


To  facilitate  integration  of  (7),  we  assume  vc  =  constant  along  the  streamline  of 
the  core  to  get 


cos  6  =  cos 


9o  + 


8 

2v. 


y2- 


y  sin  e 


de 

By 


y=0 


(8) 


This  neglects  the  point  that  vc  ^  constant  along  a  streamline,  if  there  are 
variations  in  radius  of  curvature.  Blake  (1980)  has  demonstrated  the  validity 
of  (8)  under  less  restrictive  assurptions . 

We  can  estimate  the  north-south  diameter  of  an  eddy  measured  between 
speed  maxima,  if  we  take  0O  =  tt  at  the  southernmost  extent  (the  origin)  and  0  =  J0 
at  the  northernmost  extent.  Shortly  before  an  eddy  breaks  of f  fran  the  loop  Current,  the 
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flow  is  approximately  westward  at  the  inflection  point  between  the  cyclonic  flow 
near  the  inflow  port  and  the  anti-cyclonic  flow  of  the  forming  eddy.  For  an 
infinite  radius  of  curvature  at  the  origin  (8)  becomes 

d  =  2r  =  2(vc/S)^  .  (9) 

where  d  is  the  desired  diameter  and  r  is  the  radius.  This  is  consistent  with  a 
scale  for  beta  fcossby  number,  Rg  =  vc/(0r2)  =  1.  The  penetration  of  the  current 
into  the  basin  can  be  estimated  from  (8)  by  setting  0  =  0  at  the  northernmost 
extent.  This  yields 

2  r  4  i1/2 

b  =  y-  +  4  +  2r2  (1  -  cos  0  )  I  ,  (10) 

0  i 

L°  j 

where  b  is  the  northernmost  penetration  of  the  core  of  the  current  measured  from 
the  point  of  origin,  and  rQ  is  the  radius  of  curvature  of  the  current  at  the  point 
of  the  origin.  The  numerical  model  verifies  the  usefulness  of  the  results  from 
this  crude  analysis. 

We  tested  the  validity  of  (9)  as  an  estimate  of  the  north-south  diameter  of 
the  eddies  measured  between  the  speed  maxima.  Averaged  over  35  cases,  (9) 
overestimated  the  diameter  by  6%,  when  the  maximum  speed  at  inflow,  v^n,  was  used 
to  estimate  vc.  The  eddy  diameter,  <5,^ ,  was  measured  as  the  eddy  center  passed  a 
longitude  110  km  west  of  the  western  boundary  of  the  inflow  port.  In  most  cases 
the  error  was  <5%.  In  the  cases  with  low  Reynolds  number  (see  Table  2)  or  at  very 
low  latitudes,  (9)  overestimated  the  diameter  by  more  than  10%.  The  linear 
correlation  between  the  eddy  diameters  predicted  by  (9)  and  by  the  numerical  model 
was  0.87. 

We  did  not  attempt  a  direct  verification  of  (10) .  Instead,  we  performed  a 
simpler  test  using  a  regression  equation  of  the  form 

b  =  Aq  +  A1  d  ,  (11) 
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where  Aq  and  A^  are  regression  coefficients.  This  is  an  appropriate  regression 
equation  for  (10) ,  if  all  cases  use  the  same  origin,  and  at  the  origin  0O  is  the 

same  for  every  case,  and  rQ  =  » .  Omitting  the  case  with  inflow  27°  west  of 
normal,  this  regression  equation  was  tested  using  values  of  dm  and  frcm  34  numerical 
solutions.  The  linear  correlation  between  the  penetration  distance  and  the  eddy 
diameter  was  0.986,  and  the  average  error  of  the  regression  fit  was  2%.  The  values 
tor  the  regression  coefficients  were  Aq  =  108.4  km  and  A}  =  1.232.  Note  only  1/5 
of  the  total  penetration  is  given  by  ^  and  A^>1.  When  (9)  and  the  maximum  speed 
at  inflow,  vin,  were  used  to  estimate  d  for  use  in  (11) ,  the  linear  correlation  was 
0.89,  and  the  average  error  of  fit  was  5%.  The  values  for  the  regression 
coefficients  became  Aq  =  88.6  km  and  Aj  =  1.211.  This  analysis  clearly 
demonstrates  the  importance  of  differential  rotation  in  determining  the  diameter  of 
the  eddies,  and  the  penetration  of  the  model  Loop  Current  into  the  Gulf. 

One  of  the  most  intriguing  aspects  of  the  eddy  shedding  by  the  Loop  Current  is 
the  prediction  of  the  eddy  shedding  period.  We  have  found  that  when  parameters 
realistic  for  the  Gulf  of  Mexico  are  used,  the  eddy  shedding  period  predicted  by 
the  numerical  model  is  close  to  the  quasi-annual  period  observed,  even  though  no 
time  variations  in  the  prescribed  inflow  profile  for  ^  were  included.  Thus,  we 
must  ascribe  the  eddy  shedding  to  some  instability  inherent  in  the  internal 
dynamics  of  the  system,  and  we  might  describe  the  eddy  shedding  period  as  the 
period  of  the  instability.  However,  the  Loop  Current  does  not  necessarily  meet  any 
instability  criterion  during  most  of  the  period  between  the  shedding  of  eddies. 

That  is,  the  eddy  shedding  period  is  not  determined  by  the  growth  rate  of  the 
horizontal  shear  instability,  but  by  the  time  required  for  the  Loop  Current  to  move 
into  an  unstable  configuration,  as  it  first  penetrates  into  the  Gulf  and  then  bends 
westward.  When  it  finally  bends  into  an  unstable  configuration  and  an  eddy  is 
about  to  break  off,  contours  of  growth  rate  for  the  instability  are  crossed  so 
rapidly  that  the  growth  of  the  instability  itself  contributes  little  to  the  total 
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eddy  shedding  period.  This  is  a  picture  of  a  spin-up  process  which  eventually 
becomes  unstable.  The  spin-up  is  reinitiated  after  the  break-off  of  each  eddy. 
This  spin-up  is  not  inherently  unstable  when  dissipation  is  included,  as  will  be 
shown  in  the  next  sub-section.  In  some  regions  of  parameter  space  the  northward 
penetration  and  westward  bending  proceed  to  a  steady  state,  and  no  eddy  shedding 
occurs. 

The  northward  penetration  and  westward  bending  can  be  illustrated  by 
considering  internal  Rossby  wave  dynamics  and  Reid' s  argument  in  the  form  of  a 
quasi- isostatic  adjustment  process.  For  geostrophic  divergence,  the  continuity 
equation,  (2),  for  the  reduced  gravity  model  becomes 


9h  _  0V  _  3  h  .  _3Jt_  _ 

at  T~  ~  3 1  inr  3X 


(12) 


where  Cinr  is  the  non-dispersive  internal  Rossby  wave  speed.  This  illustrates  the 
role  of  internal  Rossby  wave  dynamics  in  the  westward  spreading  of  the  Loop 
Current.  We  have  already  demonstrated  the  high,  0.991,  linear  correlation  between 
the  westward  movement  of  free  eddies  in  the  numerical  model  and  the  non-dispersive 
internal  Rossby  wave  speed.  The  maximum  rate  of  pycnocline  deepening  on  the  west 
side  of  the  Loop  implied  by  (12)  is  about  5  nv/day  for  our  standard  reduced  gravity 
case.  Now  pick  the  first  inflection  point  on  the  stream  after  inflow  on  each  of 
the  panels  in  Figure  4.  Measured  frcm  this  point,  b  from  (10)  reduces  to 


b  =  (2  ^  (1  -  cos  eQ)  )  1/2  .  (13) 

In  Figure  4  we  note  the  inflection  point  moves  only  slightly  northward  during  the 
northward  penetration  and  subsequent  westward  bending  of  the  Loop  Current.  During 
this  process  the  current  at  the  inflection  point  rotates  counterclockwise  due  to 
the  influence  of  internal  Rossby  wave  dynamics  on  the  loop,  starting  frcm  a  nearly 
northward  direction,  60=tt/2 .  As  0O  increases  frcm  tt/2  to  it,  b  increases  most  rapidly 
at  first.  Later  the  westward  bending  increases  more  rapidly.  Hence,  purely  frcm  a 
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consideration  of  Reid' s  argument  and  the  influence  of  internal  Rossby  wave 
dynamics,  we  have  a  picture  of  the  Loop  Current  first  primarily  penetrating  into 
the  basin  and  later  primarily  bending  westward.  This  is  not  an  inherently  unstable 
evolution  as  we  will  illustrate  belcw,  and  again  in  the  next  subsection.  Also, 
without  the  Rossby  wave  influence  the  current  can  assume  a  steady  source-sink 
configuration  with  no  westward  spreading  such  as  shown  in  Figure  6. 

Figure  7  shows  what  happens  when  8  is  turned  off  just  before  the  Loop  Current 
is  about  to  shed  an  eddy.  In  Figure  7  8  was  reduced  to  zero  with  a  time  constant 
of  10  days  starting  at  day  1300  for  the  standard  reduced  gravity  case  depicted  in 
Figure  5c.  Figure  7  shows  the  westward  spreading  is  halted  and  the  eddy  shedding 
process  appears  frozen.  When  8  is  reduced  to  zero  starting  20  days  later  at  day 
1320  (results  not  shown) ,  the  eddy  shedding  is  completed  as  the  source-sink  portion 
of  the  current  retreats  eastward  toward  the  steady  state  configuration.  The  eddy 
which  separates  remains  almost  stationary.  Even  apart  from  the  westward  spreading, 
8  is  not  neutral  in  the  instability  process,  since  it  influences  the  curvature  of 
the  flow.  In  baro tropic  instability  theory  for  a  mean  zonal  flow  8  is  ccrmonly 
a  stabilizing  factor.  In  our  case  it  is  clearly  a  destabilizing  factor  in  regions 
v/here  it  tends  to  increase  the  magnitude  of  the  flow  curvature. 


At  this  point  we  will  test  the  hypothesis  that  the  eddy  shedding  period  is  a 
small  multiple  of  the  time  required  for  an  eddy  to  move  an  eddy  diameter.  This 
multiple  depends  on  the  inflow  angle  of  the  current.  From  the  previous  analysis 
these  are  clearly  inportant  factors  determining  the  time  required  for  the  Loop 
Current  to  bend  into  an  unstable  configuration.  We  will  test  this  hypothesis  using 
a  regression  equation  of  the  form 

1/2  j 

pe  =  Ao  +  A1  d  0  +  cos  e0)  /  ce  (14) 
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where  pe  is  the  period  of  the  eddy  shedding,  0O  is  the  angle  of  the  inflow  with  the 
positive  x-axis,  and  ce  is  the  speed  with  which  an  eddy  moves  westward  after 
separation  from  the  Loop  Current.  In  the  first  test,  values  observed  directly  in 
the  numerical  solutions,  and  c^,,  were  used  for  d  and  ce.  The  linear  correlation 
between  the  eddy  shedding  periods  observed  in  the  numerical  solutions,  pjj,,  and  the 
predictor  in  (14)  is  0.95.  The  average  fit  error  is  8%,  and  the  regression 
coefficients  are  Aq  =  45.1  days  and  Ai  =  2.336.  Only  14%  of  the  period  is  left  to 
the  constant  AQ. 

The  regression  equation  (14)  was  tested  a  second  time  using  theoretical  values 
tor  d  and  ce.  The  values  of  d  were  obtained  from  (9)  using  Vin  to  estimate  vc. 

The  non-dispersive  internal  Rossby  wave  speed,  Cjnr,  was  used  for  ce.  The  linear 
correlation  between  predictor  and  predictand  is  0.96  with  an  average  fit  error  of 
7%.  The  regression  coefficients  are  Aq  *  56.8  days  and  A^  =  2.906.  Only  1/6  of 
the  eddy  shedding  period  is  predicted  by  the  constant,  Aq.  Of  the  35  cases,  26  had 
fit  errors  less  than  the  mean  percentage  error.  The  low  viscosity  case  with 
A=3xl06  cm2 sec" 1 ,  the  30  Sv  case  with  standard  basin  size,  and  two  of  the  narrow, 

80  km,  inflow  port  width  cases  exhibited  by  far  the  largest  fit  errors,  between  62 
and  72  days.  The  predicted  period  was  too  short  for  the  two  former,  too  long  for 
the  two  latter. 

Examination  of  Table  2  shows  that  the  role  of  friction  should  be  considered. 

The  striking  result  is  that  as  the  Reynolds  number  is  increased,  so  is  the 
eddy  shedding  period.  Four  cases  show  the  period  increasing  inversely  with  the 
eddy  viscosity,  and  six  show  the  period  increasing  directly  with  the  inf  lew  port  width. 
Some  insight  into  this  perhaps  counter-intuitive  result  can  be  gained  by 
considering  some  simple  frictional  modifications  to  the  previous  inviscid 
arguments.  We  used  the  maximum  velocity  at  inflow  when  using  (9)  to  estimate  d. 
Downstream  the  speed  at  the  core  of  the  current  may  decrease  due  to  diffusion,  or 
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entrainment  may  occur.  In  the  higher  Reynolds  number  cases  important  secondary 
circulations  occur.  These  are  not  necessarily  in  phase  with  the  eddy  shedding,  and 
may  account  for  the  irregular  eddy  shedding  period  that  occurred  when  we  set 
A  *  3  x  106cro2sec“l .  As  modifications  to  the  previous  inviscid  analysis,  these 
ideas  are  consistent  with  frictional  modifications  of  the  eddy  shedding  period 
observed  in  Table  2. 

He  have  performed  regression  analyses  using  various  measures  of  frictional 
effects  as  predictors.  For  example,  a  diffusion  time  scale  for  Laplacian  friction 
can  be  written  as  L2/A.  We  used  the  inflow  port  width  for  the  length  scale,  L. 

This  was  first  used  in  a  regression  equation  of  the  form 

Pe  =  Aq  +  A]  L2/A  (15) 

When  all  35  cases  were  used,  the  linear  correlation  between  pe  and  L2/A  was  only 
0.13,  and  the  average  fit  error  was  24%.  When  the  data  was  stratified  to  remove 
cases  where  B»  fo»  vin»  or  inflow  angle  differed  significantly  from  the  standard 
case,  the  correlation  for  the  21  remaining  cases  was  0.82.  The  average  fit  error 
was  3%.  Hcwever,  the  range  of  periods  was  very  limited,  only  90  days. 

Finally,  a  test  was  conducted  using  a  multiple  regression  equation  of  the  form 
pe  =  Aq  +  A1  d  (1  +  cos  QQ)H/ce  +  A?  L2/A  (16) 

Although  fricitonal  effects  could  be  included  in  the  multiple  regression  analysis 
in  a  number  of  ways,  we  will  discuss  only  this  simple  one.  When  d  is  calculated 
from  (9)  using  Vin  for  vc  and  Cinr  is  used  for  ce,  the  partial  correlation  for  the 
first  predictor  is  0.962,  for  the  second,  0.129.  The  multiple  correlation  is  0.978 
and  the  average  fit  error  is  6%.  The  regression  coefficients  are  Aq  *  17.2  days, 
a  »  2.932,  and  A2  “  0.121.  Only  5%  of  the  eddy  shedding  period  is  contained  in 


the  constant  Aq. 


A  wide  variety  of  cases  have  been  used  in  the  regression  analyses.  Some  of 
the  changes  should  not  affect  the  eddy  diameter  or  the  eddy  shedding  period, 
according  to  the  hypotheses  we  have  tested.  We  repeated  the  standard  case  with 
double  grid  resolution  (ax  *  10  km.  Ay  =  9.375  km  for  each  dependent  variable). 

Table  2,  Case  RG27  shows  the  solution  was  little  affected.  Changing  the 
formulation  for  friction  frcm  AV^V  to  AhV^v  in  the  standard  reduced  gravity  case 

had  only  minor  effects,  the  most  notable  being  a  10%  reduction  in  the  eddy 
diameter,  dn,.  Adding  the  Florida  Straits  to  the  model  domain  had  virtually  no 
effect,  vindicating  our  usual  treatment  of  the  eastern  port.  Changing  from  no  slip 
to  free  slip  boundary  conditions  also  had  little  effect  on  the  eddies  or  the  eddy 
shedding  except  in  the  northwest  portion  of  the  basin,  where  the  eddies  in  the  no 
slip  case  experienced  their  final  decay.  In  the  free  slip  case  the  eddies  decayed 
more  slowly  in  that  region  and  formed  northward  and  eastward  boundary  currents  in 
the  northwest  part  of  the  basin.  In  another  case  we  changed  the  shape  of  the 
prescribed  inflow  profile  for  V  from  our  usual  parabolic  shape  to  uniform  across 
the  port.  This  case  fit  well  with  the  hypotheses  we  have  tested,  suggesting  that 
the  maximum  velocity  at  inflow  is  of  primary  importance,  not  the  shape  nor  vorticity 
distribution  of  the  inflow  profile. 

By  far  the  largest  fit  error  in  the  multiple  regression  analysis,  (16) ,  was 
the  73  day  under-prediction  for  the  30  Sv  inflow  case  with  the  standard  sized  basin 
(Table  2,  Case  RG21)  when  Vin  was  used  for  vc  in  (9)  and  Cinr  for  ce.  Extending 
the  basin  300  km  northward  reduced  the  under-prediction  to  50  days,  the  second 
largest  error.  Extending  the  basin  300  km  northward  and  160  km  eastward  reduced 
the  fit  error  to  5  days.  Such  extensions  to  the  basin  had  only  minor  effect  in  the 
standard  20  Sv  case  (see  Table  2,  Cases  RG24-26).  Similar  results  were  obtained 
when  djj,  was  used  for  d  and  was  used  for  ce  in  (14) ,  shewing  the  fit  errors  are  not 


due  to  the  use  of  theoretical  values  in  (16) .  The  case  with  6  *  10"13an-lsec-l 
exhibited  even  larger  eddies  than  the  30  Sv  cases.  However,  they  were  less 
vigorous  and  their  properties  fit  well  in  most  of  the  regression  analyses.  The 
effects  of  different  port  configurations  are  discussed  in  Section  4. 

The  view  of  the  instability  of  the  eddy  shedding  Loop  Current,  presented 
earlier  in  this  section,  is  not  the  only  useful  one.  Figure  8  a,  b  offers  a 
different  picture  of  the  eddy  evolution  and  propagation.  Figure  8a  shows  the  mean 
over  3  eddy  cycles  for  the  standard  reduced  gravity  case.  This  shows  that  in  the 
mean  the  model  eddies  produce  a  western  boundary  current  in  the  Gulf  with  northward 
flow  in  the  upper  ocean.  Figure  8b  is  a  snapshot  of  the  deviation  from  the  mean  at 
day  1800.  The  deviations  depict  a  wavetrain  originating  northeast  of  the  inflow 
port  and  decaying  near  the  western  boundary.  Note  the  properties  of  the  positive 
and  negative  deviations  are  quite  similar.  This  suggests  that  for  our  nearly 
circular  eddies,  the  appropriate  regression  coefficients  for  (14)  are  Aq  =  0  and 
A1  *  2,  if  we  use  the  east-west  diameter  of  the  deviations  of  a  particular  sign  for  d 
and  Cjn  for  ce.  This  is  roughly  in  accord  with  the  values  AQ  =  45  days  and  =  2.34 
we  obtained  using  for  d.  The  east-west  diameter  of  the  deviations  for  the  standard 
reduced  gravity  case  is  d^  =  420  km  for  the  downward  (positive)  deviations 
associated  with  the  eddy  and  d2  -  505  km  for  the  upward  (negative)  deviations 
between  the  eddies,  compared  to  373  km  for  <%  from  Table  2,  Case  RGl .  To  check  the 
consistency  of  our  data,  the  eddy  shedding  period  implied  by  p>e  =  (d1+d2)/cn,  is  334  days, 
compared  to  the  327  days  actually  measured  (Table  2,  Case  RGl) . 

c.  Stabilization  of  the  Loop  Current  by  reducing  the  Reynolds  number 

We  have  already  demonstrated  that  setting  6*0  has  a  powerful  stabilizing 
effect  on  the  Loop  Current.  In  section  4  we  discuss  the  important  effects  of  port 
separation  and  topography  on  the  stability  of  the  Loop  Current.  In  this  subsection 
we  demonstrate  a  stable  regime  for  the  Loop  Current  obtained  by  reducing  the 
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Reynolds  nunber,  Re  =  where  we  take  vc  as  the  speed  at  the  core  of  the 

current,  and  as  half  the  port  width.  We  reduced  Re  in  two  ways  from  the 
standard  reduced  gravity  case,  1)  by  increasing  the  eddy  viscosity  by  a  factor  of  3 
to  A  *  3  x  107an2sec“l ,  and  2)  by  reducing  the  inflow  transport  by  a  factor  of  200 
to  obtain  an  essentially  linear  solution.  Since  the  reduced  gravity  model  required 
only  lh  minutes  of  computer  time  per  year  when  At  =  3  hours,  reducing  the  transport 
was  a  more  expedient  method  of  obtaining  a  linear  solution  comparable  to  the 
standard  case  than  linearizing  the  numerical  model.  In  the  linear  case  the 
amplitude  of  the  inflow  affects  only  the  amplitude  of  the  solution.  Linearity  was 
verified  by  comparing  solutions  with  inflows  of  1  Sv  and  0.1  Sv. 

The  lower  panels  of  Figure  8  shew  the  pycnocline  anomaly  for  (c)  the  steady 
state  solution  obtained  with  A  =  3  x  107an2sec“^>  and  (d)  the  steady  state  linear 
solution  corresponding  to  the  standard  case,  obtained  by  using  an  inflow  trans¬ 
port  of  0.1  Sv.  The  westward  spreading  due  to  differential  rotation  is  clearly 
evident  in  both  cases.  Although  the  westward  spreading  is  a  strong  destabilizing 
influence,  it  does  not  necessarily  lead  to  an  unstable  configuration,  when 
dissipation  is  included.  When  the  simple  Laplacian  friction  is  just  great  enough 
to  prevent  the  large  amplitude  discrete  eddies  we  obtained  in  the  standard  case, 
it  demonstrates  remarkable  skill  in  parameterizing  the  mean  over  an  eddy  cycle 
(compare  Fig.  8a  and  8c).  In  another  numerical  experiment  we  initialized  from 
the  steady  state  shown  in  Figure  8c.  We  kept  A  =  3  x  10  an  sec  ,  but  reduced 
0  to  zero  with  a  time  constant  of  10  days.  When  we  did  this,  the  Loop  Current 
retreated  toward  a  source-sink  flow  similar  to  Figure  6.  In  the  process  a  very 
weak  eddy  was  left  behind. 


d.  Barotropic  flat  bottom  model  results 

The  barotropic  and  reduced  gravity  models  are  mathematically  identical  but 


have  differing  gravitational  accelerations.  Hie  vorticity  equation  for  both  models 
is 

Ct  +  NL  +  f(ux  +  Vy)  +  8  v  *  A  V2c  (17) 

where  c,  is  the  relative  vorticity  and  NL  represents  the  nonlinear  terms.  In  the 
case  of  a  flat  bottom  Gulf  of  Mexico  model,  both  have  the  same  steady  state  linear 
solution  for  the  same  inflow  velocity,  except  that  the  amplitude  of  the  barotropic 
and  reduced  gravity  "h"-ananaly  fields  is  in  the  ratio  of  g'/g-  Hie  vorticity  balance 
is  between  planetary  vorticity  advection  (8v),  and  horizontal  diffusion  (AV  rJ .  For 
the  parameters  of  the  standard  case  given  in  Table  1,  the  flat  bottom  barotropic 
model  evolves  to  a  steady  state  with  a  pattern  of  free  surface  anomaly  contours 
almost  identical  to  the  linear  solution  shown  in  Figure  8d,  cf.  Figure  23a.  Some 
important  effects  of  topography  in  this  case  are  presented  in  Section  4. 

Hie  flat  bottom  barotropic  model  begins  to  exhibit  eddy  shedding  at  Reynolds  numbers 
slightly  higher  than  in  the  reduced  gravity  case.  Figure  9a  shows  the  steady  state 
solution  for  the  barotropic  model  with  a  Reynolds  number  comparable  to  a  reduced 
gravity  case  which  was  weakly  unstable.  Increasing  the  Reynolds  number  by  1/3 
produces  the  weak  eddy  shedding  regime  shown  in  Figure  9b,  with  weak  eddies 
superimposed  on  a  mean  ridge.  The  eddies  propagate  westward  along  the  ridge.  The 
reduced  gravity  model  and  the  two-layer  model  (both  with  a  flat  bottom  and  the 
topography  of  Figure  3)  exhibit  similar  regimes  with  a  weak  horizontal  shear 
instability.  However,  when  we  increased  the  Reynolds  number  further,  the  character 
of  the  barotropic  solutions  was  somewhat  different  from  the  other  models.  The 
other  models  developed  discrete  eddy  shedding  with  weak  secondary  circulations  (see 
Figures  4  and  5) .  Only  when  the  Reynolds  number  was  increased  still  further  did 
the  secondary  circulations  become  important.  In  the  barotropic  model  we  did  not 
find  the  discrete  eddy  shedding  regime  with  weak  secondary  circulations,  instead 
the  secondary  circulations  became  important  at  Reynolds  numbers  lower  than  those  at 


which  the  discrete  eddy  shedding  occurred  (see  Figs.  9c  and  9d) . 

A  shorter  baro tropic  eddy  shedding  period  (Table  2)  and  the  preceding  discuss: or. 

demonstrate  significant  differences  between  the  barotropic  flat  bottom  and  reduced 

gravity  models  due  to  the  role  of  horizontal  divergence.  Its  role  is  uninportant  in  th 

former,  important  in  the  latter.  In  the  reduced  gravity  model  the  radius  of  the  eddies 

1/2 

is  several  times  greater  than  the  internal  radius  of  deformation,  (g'h^)  /f-45  km;  in 

the  barotropic  case  several  times  less  than  the  barotropic  radius  of  deformation, 

1/2 

(gh)  /l  -  1000  km  for  h  =  300  m  as  in  Figure  9b,c,d.  Thus,  unlike  the  reduced 
gravity  model  the  horizontal  divergence  in  the  barotropic  model  is  not 
approximately  geostrophic,  and  the  Rossby  wave  propagation  is  dispersive,  as 
opposed  to  nearly  non-dispersive  in  the  reduced  gravity  case,  cf .  (12) .  In  the 
barotropic  model  a  rigid  lid  would  eliminate  the  horizontal  divergence  entirely  in 
the  flat  bottom  case,  yielding  an  infinite  barotropic  radius  of  deformation  and 
purely  dispersive  Rossby  waves. 

The  interference  of  the  secondary  circulations  and  the  distortion  of  the 
eddies  prevent  us  from  making  satisfactory  measurements  of  c^,  etc.  in  the 
barotropic  model.  This  hinders  quantitative  comparison  with  the  reduced  gravity 
model.  However,  a  qualitative  comparison  is  possible.  When  we  setB=0,  the 
solution  is  stabilized  as  in  the  other  models,  and  a  source-sink  flow  similar  to 
Figure  6  develops.  Further,  as  predicted  by  equations  (9)  and  (10)  ,  the  diameter 
of  the  eddies  and  their  penetration  into  the  Gulf  is  similar  (see  Table  2) .  Since 
the  barotropic  Rossby  waves  are  dispersive,  the  barotropic  spin-up  is  more  rapid  and 
the  eddy  shedding  period  is  shorter,  about  50  days  for  the  cases  shown  in  Figures 
9c  and  9d.  The  complicated  patterns  shown  in  these  figures  are  repeated  in  detail 
on  this  50-day  cycle.  The  stability  regimes  for  the  Loop  Current  are  re-examined 
in  section  4,  when  the  role  of  port  separation  is  included. 
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e.  Cases  of  baroclinic  or  mixed  stability 

The  standard  two-layer  flat  bottom  case  (Table  1)  exhibited  regular  discrete 
eddy  shedding  much  like  Figure  4  with  little  sign  of  other  oscillations  (Fig.  5b). 
When  the  eddy  viscosity  for  this  case  is  reduced  by  a  factor  of  3  to 
A  *  3  x  10^an^aec~l ,  the  character  of  the  eddy  shedding  is  markedly  altered. 

Figure  10  shows  a  sequence  of  six  snapshots  covering  only  one  quarter  of  an  eddy 
shedding  cycle.  A  model  representation  of  the  Florida  Straits  has  been  added  to 
the  basin.  This  did  not  significantly  affect  the  results  as  discussed  below. 

Figure  10a  looks  much  like  the  standard  case,  but  the  rest  of  the  sequence  is  quite 
different.  Figure  10b  sh<x#s  the  eddy,  which  was  about  to  shed  in  Figure  10a, 
becoming  severely  distorted.  In  Figure  10c, d  it  has  fractured  into  two  eddies. 
Figure  10e  shows  the  two  eddies  coalescing.  In  Figure  10f  they  have  returned  to 
the  form  of  a  single  eddy.  During  the  3/4ths  of  the  eddy  cycle  not  shown,  the  Loop 
Current  penetrates  into  the  Gulf,  bends  westward,  and  begins  to  form  an  eddy 
structure.  These  phases  are  much  like  the  standard  case,  except  for  a  strong 
undulation  of  the  Loop  Current. 

Figure  11a  shows  the  kinetic  energy  for  each  layer  in  the  rectangular  part  of 
the  basin  as  a  function  of  time.  Clearly  evident  are  two  fairly  regular 
oscillations,  one  of  380  days,  the  other  of  53  days.  The  longer  one  is  the 
eddy  shedding  cycle.  The  shorter  is  associated  with  an  undulation  of  the  Loop 
Current  and  the  fracturing  of  the  eddies.  Because  these  two  cycles  are  not 
harmonically  related,  the  details  of  each  eddy  cycle  are  quite  different.  Figures 
lib  and  11c  show  the  absence  of  the  high  frequency  oscillation  in  the 
standard  two-layer  flat  bottom  case,  and  in  the  reduced  gravity  case  with 
A  *  3  x  106cm2sec-1  (Table  2,  Case  RG8) .  The  latter  is  the  same  eddy  viscosity 
used  for  the  case  shown  in  Figures  10  and  11a.  Another  difference  in  Figure  11a 
from  the  standard  case  (Fig.  lib)  is  the  periodic  large  growth  of  the  kinetic 
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energy  in  the  lc*rer  layer,  lagging  the  peaks  in  the  upper  layer  by  about  90°.  The  time 
that  the  two  energy  curves  in  Figure  11a  coincide  corresponds  to  Figures  10d  and  10e, 
showing  the  baroclinic  to  barotropic  energy  conversion  which  takes  place  as  the 
eddies  fracture.  Based  on  the  results  of  McWilliams  and  Flierl  (1979) ,  we 
anticipate  that  Rossby  wave  dispersion  in  the  lower  layer  plays  an  important 
catalytic  role  in  the  subsequent  rapid  decrease  in  the  lower  layer  kinetic 
energy. 

Twice  we  have  repeated  the  experiment  with  the  energy  curves  shown  in  Figure 
11a  without  including  the  Florida  Straits  as  an  extension  to  the  basin,  once  with 
and  once  without  bottom  friction.  In  both  cases  the  energy  curves  matched  those 
shown  in  Figure  11a  in  detail  for  over  two  years.  After  that  they  diverged 
quantitatively,  but  not  qualitatively.  The  energy  curves  for  the  cases  with  and 
without  bottom  friction  agreed  in  detail  for  the  three  year  duration  of  the 
comparison.  The  bottom  friction  was  quadratic  with  a  drag  coefficient  of 
1.3  x  10”3.  The  other  flat  bottom  experiments  we  performed  should  be  even  less 
sensitive  to  bottom  friction.  We  repeated  the  experiments  discussed  above,  but 
included  the  topography  of  Figure  3  and  increased  the  inflow  in  the  upper  layer  to 
25  Sv.  Seme  short  term  fluctuations  (~50  days)  remained  in  the  energy  curves  when 
there  was  no  bottom  friction.  These  were  reduced  or  eliminated  when  the  bottom 
friction  was  added,  although  they  were  also  removed  by  increasing  the  lateral 
fr iction. 

Comparison  with  the  reduced  gravity  results  for  the  same  eddy  viscosity  (Fig. 

11c,  Table  2,  Case  RG8)  suggests  that  the  two-layer  case  (Figs.  10,  11a)  is 
exhibiting  a  mixed  instability,  although  it  is  not  clear  in  this  case  that  a  con¬ 
tribution  from  the  horizontal  shear  instability  is  necessary  to  explain  the  eddy 
shedding.  As  discussed  in  section  3b,  it  is  not  needed  to  explain  the  Loop  Current 
penetration,  westward  spreading,  or  tendency  to  curve  back  on  itself.  The  rapid  rise  in 
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the  lower  layer  kinetic  energy  lagging  the  peak  in  the  upper  layer  kinetic  energy  is 
consistent  with  results  obtained  by  Holland  and  Lin  (1975)  for  mid-latitude  mesoscale 
eddies.  The  period  of  the  high  frequency  oscillation  is  also  similar. 

To  help  clarify  the  analysis,  we  looked  for  a  case  with  a  pure  signature  of 
baroclinic  instability.  Figure  12  shows  results  for  a  two-layer  flat  bottom  case 
with  the  inflow  angled  27°  west  of  normal  and  confined  to  the  upper  layer.  Hie 
snapshots  of  pycnocline  anomaly  30  days  apart  (Fig.  12a, b)  look  much  like  the  weak 
horizontal  shear  instability  regime.  However,  the  phase  plot  of  pycnocline  anomaly 
versus  time  at  x  =  1010  km  from  the  western  boundary  (Fig.  12c)  shews  the  period  is  much 
shorter,  58  days.  Furthermore,  there  is  little  evidence  of  modulation  by  a  longer 
time  scale,  except  for  the  first  eddy.  The  amplitude  of  the  oscillation  is  much 
stronger  on  the  westward  flowing  south  side  of  the  loop,  consistent  with  a 
baroclinic  instability;  e.g..  Philander  (1976),  Rhines  (1977). 

Figure  13a  shows  the  kinetic  energy  in  the  basin  as  a  function  of  time  in  the 
upper  and  lower  layers.  It  is  very  similar  to  the  results  obtained  for  baroclinic 
instability  by  Holland  and  Lin  (1975) ,  including  the  maximum  in  upper  layer  kinetic 
energy  near  the  onset  of  baroclinic  instability  followed  by  a  rise  in  lower  layer 
kinetic  energy.  The  period  of  the  oscillation  is  also  similar.  In  contrast  the 
corresponding  reduced  gravity  case  exhibits  discrete  eddy  shedding  similar  to 
Figure  4  with  a  period  of  284  days  (Table  2,  Case  RG39) .  For  comparison  the 
kinetic  energy  in  the  basin  for  the  reduced  gravity  case  is  plotted  as  a  function 
of  time  in  Figure  13b.  A  detailed  investigation  of  eddy  energetics  in  the  Gulf  of 
Mexico  is  in  preparation. 

Figure  14  shews  spectra  of  the  pycnocline  ancmaly  at  a  particular  location  in  the 
basin  for  four  different  cases,  a)  the  standard  two-layer  flat  bottem  case,  b)  the 
reduced  gravity  case  with  A  =  3  x  lB^cm^sec-!,  c)  the  two-layer  flat  bottom  case 

with  A  =  3  x  10^an2sec~l,  and  d)  the  two- layer  flat  bottem  case  with  inflow  27° 
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west  of  normal  and  confined  to  the  upper  layer.  Parameters  not  mentioned  have  the 
same  values  as  the  standard  case  (Table  1) .  Figures  14a, b,c  exhibit  a  major  peak 
in  the  vicinity  of  one  year  corresponding  to  the  discrete  eddy  shedding  cycle  found 
in  the  standard  case.  In  Figures  14a, b  no  other  major  peaks  are  found,  except  for 
harmonics.  This  occurs  despite  the  irregular  shedding  period  found  in  the  case 
shown  in  Figure  14b.  The  irregular  shedding  is  exhibited  by  higher  energy  levels 
at  higher  frequencies  and  a  noisier  looking  spectrum,  c.  f.  Fig.  30a.  In  contrast 
Figure  14d  exhibits  a  major  peak  only  at  6  cycles/year.  A  peak  near  6  cycles/year 
is  also  found  in  Figure  14c. 

Earlier  we  noted  that  the  discrete  eddy  shedding  regime  associated  with  the 
horizontal  shear  instability  most  resembles  observations.  Since  baroclinic 
instability  becomes  important  when  we  reduce  the  eddy  viscosity  in  the  two-layer 
flat  bottom  model,  does  this  negate  the  value  of  the  reduced  gravity  results?  When 
baroclinic  instability  becomes  important,  the  results  appear  less  realistic.  This 
question  is  re-examined  in  section  4  when  we  include  the  topography  shown  in 
Figure  3. 
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4.  External  influences  on  the  eddy  shedding  by  the  Loop  Current 


We  have  shown  that  several  distinct  flow  regimes  can  exist  in  the  models' 
statistically  steady  state.  While  experiments  dominated  by  a  horizontal  shear 
instability  of  the  internal  mode  produced  the  eddy  shedding  behavior  most  closely 
resembling  that  from  the  actual  Loop  Current,  baroclinic  instability,  mixed 
instabilities,  and  purely  steady  flows  have  also  been  observed  in  individual 
experiments  with  the  three  model  prototypes,  within  the  shear  instability  regime, 
variation  of  parameters  lead  to  different  periods  for  eddy  shedding.  Constant 
inflow  transport  and  realistic  values  for  the  parameters  yield  an  eddy  shedding 
period  near  annual.  This  suggests  that  in  the  actual  Loop  Current,  regular, 
quasi-annual,  eddy  shedding  might  occur  with  no  significant  variations  in  inflow 
transport  or  velocity  profile  —  the  shedding  being  controlled  by  internal  dynamics 
within  the  Gulf.  To  further  examine  this  hypothesis  it  is  necessary  to  test  the 
sensitivity  of  the  models  to  spatial  or  temporal  changes  in  the  boundary  conditions 
and  to  the  inclusion  of  realistic  topography. 

In  this  section  we  first  describe  the  response  of  the  model  to  changes  in 
location,  width,  and  separation  of  the  inflow  and  outflow  ports.  Results  from 
these  experiments  and  those  in  section  3  are  summarized  in  a  stability  regime 
diagram  using  appropriate  definitions  of  the  Reynolds  number  and  beta  Rossby  number 
as  the  axes.  Next,  the  dynamical  influence  of  the  large  amplitude  topography  in 
the  Gulf  of  Mexico  is  analyzed.  Particular  attention  is  given  to  the  topography  of 
the  West  Florida  Shelf.  In  the  presence  of  deep  water  inflow  through  the  Yucatan 
Straits,  we  have  found  that  this  topography  exerts  a  crucial  influence  on  the  Loop 
Current  penetration  and  subsequent  eddy  shedding.  Finally,  we  explore  the  model 
response  to  spatial  and  tanporal  variations  in  the  inflow  transport. 

Most  experiments  in  this  section  utilize  the  reduced  gravity  model  since  it  is 
less  costly  than  the  two- layer  model,  but  still  reproduces  the  growth  of  the  Loop 
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Current  and  generation  of  Loop  eddies  characteristic  of  the  horizontal  shear 
instability  regime.  Experiments  involving  bottom  topography  utilize  both  the 
barotropic  and  two-layer  models. 

a.  Port  location,  boundary  conditions,  and  irregular  coastlines 

Section  3  has  clearly  demonstrated  that  a  simple  potential  vorticity 
conservation  model  of  the  Loop  Current  is  quite  useful  for  interpreting  our 
results.  However,  we  might  anticipate  that  the  flow  regime,  the  eddy  shedding 
period,  the  diameter  of  the  eddies,  the  extent  of  Loop  Current  penetration  into  the 
Gulf,  and  other  features  of  the  model  solutions  also  exhibit  some  dependence  on  the 
width  and  separation  of  the  ports. 

The  standard  20  Sv  reduced  gravity  case  was  repeated  with  the  southern  port 
moved  200  km  westward  and  the  eastern  port  150  km  northward.  Figures  5c  and  15 
compare  the  standard  experiment  and  the  moved  port  case  using  a  t-y  plot  of  PHA  190  km 
west  of  the  southern  port  in  each  case.  In  the  standard  experiment  the  first  eddy 
passes  with  maximum  amplitude  at  day  335.  Subsequent  eddies  pass  at  an  interval  of 
327  days.  In  the  moved  ports  experiment  the  first  eddy  passes  with  maximum 
amplitude  at  day  285  with  subsequent  eddies  shed  at  virtually  the  same  interval  as 
before  (Table  2,  Cases  RG1  and  RG35) . 

If  we  account  for  the  phase  difference  of  the  first  eddy  due  to  spin-up,  we 
obtain  a  direct  comparison  of  spatial  structure  for  the  two  experiments  as  shown  in 
Figure  16.  These  PHA  snapshots  taken  50  days  apart,  but  in  similar  phases  of  model 
evolution,  dramatically  illustrate  the  insensitivity  of  the  eddy  shedding  behavior 
to  changes  in  port  location  or  separation.  The  eddy  diameter,  the  latitude  of 
northernmost  penetration  of  the  Loop  Current,  and  the  eddy  shedding  rate  have  been 
affected  very  little. 

If  the  inflow  port  is  moved  still  further  to  the  west,  eddy  shedding  occurs  in 
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a  manner  qualitatively  similar  to  that  shown  in  Figure  16.  However,  the  western 
boundary  begins  to  interfere  with  the  eddy  shedding  process  and  slows  the  shedding 
rate.  For  example,  when  the  inflow  port  is  centered  280  km  from  the  western 
boundary,  the  eddy  shedding  period  increases  to  390  days.  Moving  the  port  still 
further  west  leads  to  multiyear  eddy  shedding  periods.  If  the  two  ports  are  moved 
as  close  together  as  possible,  with  the  southern  port  moved  as  far  east  as  possible 
and  the  eastern  port  moved  as  far  south  as  possible,  the  flow  sinply  enters  from 
the  south  and  exits  immediately  to  the  east.  No  Loop  Current  growth  or  eddy 
shedding  is  observed. 

In.  Fig.  17  we  effectively  reduce  the  port  separation  by  inserting  a  land  mass. 
Results  fran  these  moved  ports  experiments  demonstrate  that  if  the  ports  are  closer 
than  a  certain  mimimum  distance,  a  steady  source-sink  regime  with  no  westward  spreading 
replaces  the  eddy  shedding  regime.  A  similar  result  occurs  when  we  set  8=0  (Fig.  6). 
Even  if  the  ports  are  sufficiently  separated  to  allow  eddy  shedding,  results  in  section 
3  shew  that  reducing  the  Reynolds  number  sufficiently  will  prevent  the  eddy  shedding 
and  produce  a  steady  regime  with  westward  spreading  of  the  Loop  (Fig.  8c) . 

To  synthesize  these  diverse  model  results,  it  is  useful  to  introduce  two 
nondimens ional  parameters  based  on  the  vorticity  equation,  (17) ,  and  the  model 
results:  the  Reynolds  number.  Re  =  vinLh/A,  and  the  beta  Rossby  number, 

RB  =  vin/(BLp2) ,  where  Vj^  is  the  maximum  speed  of  the  current  on  inflow,  Lh  is 
half  the  width  of  the  inflow  port,  and  Lp  is  half  the  distance  between  the  centers 
of  the  ports.  Results  from  Table  2  for  barotropic  and  reduced  gravity 
experiments  have  been  plotted  as  a  function  of  Re  and  Rg  in  a  regime  diagram 
(Fig.  18).  The  three  regimes,  eddy  shedding  (E) ,  steady  westward  spreading  (W) , 
and  steady  source-sink  (N)  are  indicated  by  symbols  for  each  experiment.  This  is 
not  meant  to  preclude  the  existence  of  additional  regimes  at  higher  Reynolds 
numbers  where  instabilities  may  occur  on  scales  small  compared  to  the  eddies  shed 


by  the  Loop  Current  in  these  experiments.  Each  regime  is  characterized  by  a 
distinct  range  of  Re  and  Rq.  The  critical  beta  Rossby  nixnber,  Rqc  -  2  for 
transition  to  regime  N.  For  the  reduced  gravity  experiments  the  critical  Reynolds 
number,  Rec  =  25  for  transition  from  regime  W  to  regime  E.  For  the  barotropic 
cases  Rec  =  40.  Although  we  have  not  attempted  to  delineate  the  regime  boundaries 
in  detail,  several  experiments  have  exhibited  characteristics  of  two  regimes.  For 
example,  near  the  transition  from  regime  E  to  regime  W,  all  the  models  exhibited 
eddy-shedding  superimposed  on  a  westward  spreading  ridge  in  the  free  surface 
(Fig.  9b) . 

Using  (9)  from  the  conservation  of  potential  vorticity  analysis  in  section  3, 
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we  can  also  express  %  as  the  ratio  of  the  eddy  radius,  r  =  (Vin/B)  ,  to  the  port 
separation,  2  Ip.  For  approximately  r/Lp  >  v  2,  neither  westward  spreading  nor 
eddy  shedding  occur.  The  critical  port  separation  can  be  expressed  as 

Lpc  -  <*)„  /  *W1/2  ■  r  /0W,/2  <>8) 

where  R^c  ~2  and  Lpc  is  half  the  critical  distance  between  the  centers  of  the 
ports.  In  our  experiments  we  have  not  detected  any  dependence  of  Rec  on  Re  or 
other  model  parameters,  although  we  would  anticipate  some  dependence  on  the  angle 
of  inflow,  the  geometry,  and  the  orientation  of  the  basin.  The  importance  of  rb  to 
the  stability  of  the  Loop  Current  has  also  been  noted  by  T.  Sugimoto  (personal 
comnunication) ,  who  found  all  three  regimes  in  laboratory  tank  simulations  of  the 
Gulf  circulation. 

Using  the  regime  diagram  (Fig.  18) ,  the  vorticity  equation,  (17) ,  and  the 
conservation  of  potential  vorticity  analysis  in  section  3,  we  can  state  several 
general  conclusions  about  the  model  Loop  Current  and  the  eddy  shedding.  In  regime 
W  the  flow  spreads  westward  and  evolves  to  a  steady  state  with  no  eddy  shedding. 

The  flow  resembles  the  linear  viscous  solution  (Fig.  8d)  where  planetary  vorticity 


advection  ( Bv)  balances  horizontal  friction  (A^c).  weak  nonlinearity  may  be 
present,  causing  the  westward  spreading  to  occur  at  a  higher  latitude  than  if  the 
flow  were  linear  (Figs.  8c,  9a) . 


Fran  (17)  nonlinearity  becomes  important  in  regime  E.  In  the  presence  of  dif¬ 
ferential  rotation,  6^ and  increased  conservation  of  potential  vorticity,  the  Loop  Current 
tends  to  curve  back  on  itself,  cutting  off  eddies  (e.g.,  Figs.  4,  16).  *Ihe  role  of 
the  divergence  term,  f(ux  +  Vy) ,  was  discussed  in  section  3d.  Near  the  transition 
between  regimes  W  and  E  weak  eddies  are  superimposed  on  a  mean  westward  spreading 
loop  (Fig.  9b) . 

In  regime  N  where  %  =  v ^  (bl2  )  >2,  the  nonlinearity,  the  small  port 
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separation  or  small  6,  and  the  tendency  of  the  flow  to  conserve  potential  vorticity 

while  in  approximate  geostrophic  balance  all  act  in  conjunction  to  prevent  the 

westward  spreading  and  to  limit  the  northward  penetration  of  the  current.  The 

nonlinear  terms  counter  the  westward  spreading  due  to  planetary  vorticity 

advection,  the  eastern  port  fixes  the  location  of  the  current  at  outflow,  and  the 

tendency  to  conserve  potential  vorticity  influences  the  path  of  the  current.  For 

normal  inflow  with  no  curvature,  (8)  predicts  no  westward  spreading  and  (13)  predicts 
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the  maximum  northward  penetration  of  the  loop  Current  to  be  b  =  (2  Vc/B)  ,  when  it  is 
not  constrained  by  a  small  port  separation.  The  currents  in  regime  N  evolve  to  a 
steady  state  before  penetrating  that  far.  These  results  suggest  that  the  current 
must  penetrate  approximately  that  distance  for  westward  spreading  and  eddy  shedding 
to  occur.  Given  the  constraint  of  a  small  port  separation,  it  is  evident  from 
sketching  alternative  paths  that  paths  with  substantially  greater  or  lesser 
penetrations  would  conform  even  less  well  with  (8)  than  that  shown  in  Figure  17. 

In  regime  N  cases  with  low  Reynolds  numbers,  the  horizontal  diffusion  would  tend  to 
prevent  sharp  curvature  of  the  current. 

Our  experiments  thus  far  suggest  that, in  the  reduced  gravity  model  with  Re>25, 
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potential  vorticity  is  so  well  conserved  in  the  Loop  Current  that,  within  limits,  port 
location  has  little  influence  on  its  behavior.  To  explore  this  point  further,  experiments 
with  the  reduced  gravity  model  were  conducted  in  which  the  entire  eastern  boundary 
was  opened  (Table  2,  Cases  RG44-46) .  The  Loop  Current  was  free  to  choose  any 
path  out  of  the  basin  along  the  eastern  boundary.  The  open  eastern  boundary 
conditions  were  applied  as  before,  this  time  to  the  entire  eastern  side  of  the 
model.  Figure  19  provides  a  snapshot  and  a  phase  plot  of  PHA  from  an  experiment 
with  an  east-west  extent  for  the  model  domain  of  1760  km.  The  phase  plot  should  be 
compared  to  Figure  5c  and  the  snapshot  was  chosen  to  be  in  a  phase  of  eddy  shedding 
similar  to  Figure  7a.  Hie  Loop  Current  still  exits  the  basin  along  the  southern 
boundary  at  all  times.  Hie  eddy  shedding  period  is  similar  to  the  standard  20  Sv 
reduced  gravity  case  (Table  2,  Cases  RG1  and  RG44) ,  but  slowly  increases 
monotonically  as  the  Loop  Current  spins  up  a  gyre  in  the  eastern  Gulf  which  has  the 
north-south  scale  of  the  basin.  This  spin-up  occurs  on  a  time  scale  long  compared 
to  the  eddy  shedding  period. 

Based  on  these  results,  and  recognizing  the  strong  tendency  for  the  Loop 
Current  to  conserve  potential  vorticity,  we  speculate  that  the  actual  Loop  Current 
may  have  "chosen"  its  present  outflow  path  from  the  Gulf  of  Mexico.  Florida  at  one 
time  was  totally  submerged  and  the  extensive  carbonate  seabed  (Martin  and  Case, 

1975)  may  have  been  eroded  by  the  Loop  Current  over  geological  time  scales  to  form 
the  Florida  Straits. 

Though  we  have  examined  the  models'  sensitivity  to  changes  in  port  location, 
it  is  also  important  to  investigate  their  sensitivity  to  changes  in  details  of  the 
inflow  and  outflow  boundary  conditions,  in  this  case  model  response  was 
investigated  by  moving  the  specification  of  the  inflow  and  outflow  conditions  from 
inmediately  adjacent  to  the  basin  to  positions  upstream  on  inflow  and  downstream  on 
outflow.  Hui8,  the  Yucatan  and  Florida  Straits  were  added  to  the  domain  by  the 


45 


addition  of  three  rectangular  land  masses.  Figure  20  provides  results  of  this 
experiment  for  the  standard  reduced  gravity  model  and  two  additional  experiments  in 
which  the  width  of  the  Florida  Straits  was  first  halved  and  then  doubled,  while  the 
inflow  transport  remained  constant.  Figure  20  shows  the  PHA  at  day  1080  for  all 
three  experiments.  Clearly  there  is  little  qualitative  difference  in  the  eddy  size 
in  the  three  cases  and,  even  after  three  years,  very  little  phase  difference  in  the 
eddy  shedding  cycle.  (See  Table  2,  Cases  JG32-34  for  additional  canparisons) . 

Intuitively,  it  seems  likely  that  the  domain-averaged  kinetic  and  potential 
energies  would  be  considerably  different  for  the  three  different  outflow  port 
widths  considered.  For  example,  since  the  transport  on  outflow  must  equal  that  on 
inflow,  a  higher  flux  of  kinetic  energy  on  outflow  would  be  expected  for  the 
half-width  case  as  compared  to  the  other  two.  To  investigate  this  further  we 
consider  the  kinetic  energy  equations  for  the  reduced  gravity  model: 

(ke)  +  v  *  (v  ke)  +  v  •  (Up)  -  Pv-^l  -  p  v  •  f  (19) 

(PE)  +  P  v  •  V  =  0  (20) 

where  P  =  pg'(h-H),  KE  =  j  ph(u^  +  v^)  and  PE  =  pg'(h-H)^  .  Here  h  is  the 
instantaneous  layer  thickness  and  h  is  the  mean  thickness. 

The  second  term  on  the  left-hand  side  of  (19)  represents  kinetic  energy  flux 
divergence,  the  third  term  represents  pressure  work  divergence,  and  the  fourth  term 
represents  potential  to  kinetic  energy  conversion.  The  term  on  the  right  is  the 
dissipation.  Averaging  over  the  rectangular  part  of  the  domain,  contributions  to 
the  kinetic  energy  flux  and  pressure  work  occur  only  at  the  ports.  Since  we  have 
included  the  straits  we  are  able  to  calculate  flux  terms  at  the  ports  directly. 

Over  an  eddy  cycle,  in  statistical  equilibrium,  the  tendency  terms  for  the 
domain  averaged  kinetic  energy  and  potential  energy  vanish,  as  do  the  potential  to 
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kinetic  energy  conversion  terms.  Figure  21  presents  box  diagrams  for  the 
domain-averaged  energy  for  the  standard  reduced  gravity  experiment,  the  half-width 
outflow  experiment,  and  the  double  width  outflow  experiment,  respectively.  Time 

averages  of  each  term  were  taken  over  exactly  one  eddy  cycle  in  each  case. 

For  the  standard  case  the  rate  of  pressure  working  on  the  inflow  and  outflow 
is  an  order  of  magnitude  larger  than  the  kinetic  energy  flux  on  inflow  or  outflow. 
The  net  rate  of  pressure  working  on  the  basin  is  roughly  four  times  the  rate  of  net 
kinetic  energy  flux  into  the  basin.  For  the  half-width  outflow  experiment  the  rate 
of  kinetic  energy  outflow  is  quadrupled.  However,  the  pressure  work  term  on 
outflow  is  only  about  half  the  value  for  the  standard  case.  The  net  result  is  that 
the  total  kinetic  and  potential  energy  for  the  basin  remains  almost  the  same  as  in 
the  standard  case.  For  the  double-width  experiment  the  kinetic  and  potential 
energies  in  the  basin  are  only  slightly  larger  than  in  the  other  two  cases  and  the 
kinetic  energy  outflow  is  reduced  by  about  16  percent  from  the  standard  case. 

While  the  kinetic  energy  flux  and  the  rate  of  doing  pressure  work  at  the  outflow 
port  may  vary  markedly  among  these  experiments,  the  net  kinetic  and  potential 
energy  in  the  basin  varies  little. 

Before  leaving  this  aspect  of  model  response  to  external  variables  we  wish  to 
evaluate  the  influence  of  a  coastline  irregularity  in  the  actual  Gulf  of  Mexico, 
the  Yucatan  Peninsula.  To  investigate  this  external  influence  we  conducted  an 
experiment  identical  to  the  standard  reduced  gravity  case  but  with  the  Yucatan 
Peninsula  represented  as  a  rectangular  block  just  west  of  the  inflow  port.  It 
should  be  noted  in  Figure  2  that  this  geometry  actually  better  represents  the 
Yucatan  Shelf  than  the  Yucatan  Peninsula.  Figure  22  presents  a  snapshot  of  PHA  for 
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that  case  at  day  1080  and  a  phase  plot  at  x  *  1010  km.  While  the  eddy  shedding 
still  occurs,  the  eddy  shedding  period  for  this  case  is  430  days,  100  days  longer 
than  for  the  standard  case.  This  is  due  to  at  least  three  factors.  The  first,  and 
apparently  most  inportant,  is  that  while  the  minimum  penetration  latitude  does  not 
change  markedly,  the  maximum  latitude  of  Loop  Current  penetration  is  about  100  km 
greater.  This  means  that  in  the  Yucatan  Peninsula  case,  the  Loop  Current  must 
undergo  greater  changes  in  its  penetration  distance  over  an  eddy  cycle,  thus 
increasing  the  eddy  shedding  period.  In  section  3  we  noted  that  reducing  the 
north-south  extent  of  the  basin  can  in  some  cases  increase  the  eddy  shedding 
period.  This  was  effectively  done  in  part  of  the  basin  when  the  Yucatan  Peninsula 
was  inserted.  The  eddies  were  also  required  to  shed  at  a  higher  latitude,  reducing 
the  internal  Rossby  wave  speed  by  about  7%,  thus  from  (14)  potentially  increasing 
the  eddy  shedding  period  by  a  similar  amount.  However,  this  effect  was  not 
reflected  in  the  actual  movement  of  the  eddies  (Table  2,  Cases  RG1  and  HG37) . 

b.  Bottom  topographic  influences 

Thus  far  we  have  examined  the  dynamics  of  the  Loop  Current  in  terms  of  a 
reduced  gravity  system  and  a  barotropic  or  two-layer  system  with  a  flat  bottom. 
Figure  2  shews  that  the  bottom  topography  of  the  Gulf  of  Mexico  varies 
dramatically  over  the  basin.  While  the  bottom  is  deeper  than  3500  m  in  the  western 
Gulf,  it  is  180  m  deep  or  less  over  roughly  20%  of  the  Gulf  area.  Along  the 
Canpeche  escarpment  the  bottom  slopes  are  greater  than  one  percent  and  in  the 
vicinity  of  the  Yucatan  and  Florida  Straits  they  may  be  very  nearly  as  large.  From 
simple  scaling  of  the  depth-averaged  vorticity  equation,  it  is  clear  that  in  the 
vicinity  of  the  Loop  Current,  bottom  slopes  may  result  in  vorticity  advection 
(through  the  topographic  beta  effect,  -  ^  v.Vh)  larger  than  the  planetary  vorticity 
advection.  We  have  already  demonstrated  the  crucial  role  of  beta  in  the 
eddy  shedding  process.  We  must  now  investigate  the  additional  influence  of  the 


topography.  Further,  it  is  likely  that  the  baroclinic  instability  and  mixed 
instability  which  we  observed  in  the  flat  bottom  experiments  is  considerably 
modified  by  the  bottom  roughness.  Rhines  (1977)  has  shown  that  significant  bottom 
variations  on  scales  greater  than  the  deformation  radius  are  iiqportant  for 
preserving  vertical  shear.  Without  this  bottom  influence  the  energy  cascade  toward 
depth- independent  flow  may  have  a  time  scale  of  only  a  few  months.  The  results  of 
McWilliams  and  Flierl  (1979)  suggest  that  such  a  cascade  would  make  the  eddies 
increasingly  subject  to  dispersion,  particularly  in  the  lower  layer. 

We  have  included  the  Gulf  topography  in  the  idealized  fashion  discussed  in 
section  2  and  shown  in  Figure  3.  We  have  retained  the  large  scale  features  of  the 
topography,  maintained  the  proper  depth  scales  and  preserved  the  significant  slopes 
in  this  idealization.  Much  of  the  topographic  variability  on  scales  less  than  the 
deformation  radius  has  been  removed. 

To  illustrate  the  control  of  the  topography  on  the  barotropic  mode,  we 
designed  a  one-layer  experiment  with  topography  which  was  identical  to  a  flat 
bottom  barotropic  experiment  in  which  no  eddy  shedding  occurred.  The  deepest  water 
was  3000  m  in  each  case,  and  a  constant  30  Sv  inflow  was  applied  with  a  parabolic 
cross-stream  inflow  profile.  Other  parameters  are  also  standard  (Table  1) .  Figure 
23  shews  the  steady  state  FSA  for  both  cases.  No  eddy  shedding  occurred  in  either 
experiment.  Note  that  the  bottom  topography  exerts  such  overwhelming  control  on 
the  current  in  Figure  23b  that  upon  entering  the  Gulf  it  traverses  a  5000  km  path 
clockwise  around  the  basin  to  the  outflow  port  rather  than  take  the  400  km  direct 
path.  The  strong  tendency  for  this  barotropic  current  to  follow  f/h  contours  is 
reflected  in  the  outline  of  the  FSA,  which  looks  v?ry  similar  to  the  bottom  depth 
contours  in  Figure  3.  Through  the  topographic  beta  effect,  the  divergence  term  in 
the  vorticity  equation,  (17) ,  is  of  crucial  importance  in  the  case  with  topography, 
unlike  the  flat  bottom  case. 


49 


Before  investigating  direct  topographic  influences  due  to  inflow  in  the  deep 
water,  we  will  conduct  a  set  of  three  nearly  identical  experiments  with  standard 
parameters  (Table  1),  except  that  in  the  first  experiment  we  use  the  reduced 
gravity  model  with  25  Sv  inflow,  in  the  second  we  use  the  two-layer,  flat  bottom 
model  with  25  Sv  in  the  upper  layer,  no  inflow  transport  in  the  lower  layer,  and  in 
the  third  we  add  bottom  topography  to  the  second  experiment.  These  experiments 
were  integrated  to  statistical  equilibrum  for  six  years.  Figure  24  shows  snapshots 
of  PHA  for  each  experiment  during  year  6  and  Figure  25  shows  phase  plots  of  PHA  for 
the  first  five  years  in  a  north-south  section  at  x  =  810  km. 

Although  each  experiment  was  forced  with  the  same  inflow  transport  in  the 
upper  layer,  it  is  clear  that  the  energy  of  the  eddy  following  separation  is 
markedly  different  in  the  three  cases.  Loss  of  energy  to  the  deep  water  through  a 
generalized  baroclinic  instability  (Rhines,  1977)  is  impossible  in  the  reduced 
gravity  cases.  Consequently,  the  eddy  maintains  its  intensity  as  it  approaches  the 
western  boundary.  However,  in  the  flat  bottom  case,  where  loss  of  energy  to  the 
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deep  water  is  possible,  the  eddies  are  weaker  than  in  the  reduced  gravity 
experiment  and  lose  much  of  their  energy  by  the  time  they  reach  the  western 
boundary.  Further,  by  the  time  they  reach  the  western  Gulf  they  have  a  smaller 
diameter  than  the  eddies  from  the  other  two  experiments.  For  the  bottom  topography 
experiment  the  eddy  intensity  is  intermediate  between  the  reduced  gravity  and  the 
two-layer  flat  bottom  experiment.  Also,  for  both  the  topography  experiment  and  the 
reduced  gravity  experiment  the  shed  eddies  are  discrete,  whereas  in  the  flat  bottom 
experiment  we  observe  a  westward  bending  loop  with  eddies  imbedded  in  it.  The 
latter  is  characteristic  of  the  weak  eddy  shedding  regime  near  the  transition 
between  regimes  E  and  W  (Fig.  18).  The  case  with  a  pure  signature  of  baroclinic 
instability  (Figs.  12,  13a)  is  identical  to  the  flat  bottom  case  discussed  here 
except  that  the  inflow  is  27°  west  of  normal  and  the  upper  layer  inflow  transport 
is  20  Sv  rather  than  25  Sv.  It  also  exhibits  eddies  super  imposed  on  a  westward 
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bending  loop,  but  the  period  of  the  eddies  is  only  58  days. 

The  eddy  shedding  period  in  days  for  these  three  cases  is  360  for  the  reduced 
gravity,  161  for  the  two-layer  flat  bottom,  and  248  for  the  two-layer  with  Gulf 
topography.  The  short  eddy  shedding  period  in  the  two-layer  flat  bottom  case  is 
apparently  a  consequence  of  the  loss  of  energy  from  the  upper  to  the  lower  layer, 
which  tends  to  weaken  the  Loop  Current.  The  eddies  shed  from  the  Loop  are  thus 
neither  as  intense,  nor  as  large  (r= (vc/B) in  the  two- layer  flat  bottom  experiment 
as  in  the  reduced  gravity  experiment.  Since  we  found  in  section  3  that  the  eddy 
shedding  period  was,  to  a  good  approximation,  directly  proportional  to  the  eddy 
size,  (14),  we  expect  weaker  eddies  to  shed  at  a  higher  rate.  Our  results  are 
consistent  with  the  prediction  that  the  eddy  shedding  period  is  shortest  in  the 
experiment  with  greatest  energy  loss  from  the  Loop  Current  to  the  lower  layer  and 
longest  for  the  reduced  gravity  case  with  no  transfer  of  energy  to  the  lower  layer. 

Even  the  case  with  a  pure  signature  of  baroclinic  instability  is  quite 
consistent  with  these  ideas,  both  in  terms  of  the  r-  -*-/2  scale  and  an  eddy  period 

of  pe  =  2d/ce,  c.f.  (14),  where  Ag  =  0  and  A i  =  2  are  regression  coefficients 
suggested  at  the  end  of  section  3b  for  a  continuous  wave  train,  when  d  is  the 
east-west  diameter  of  an  eddy.  Thus,  in  a  mixed  stability  regime  we  find  an  eddy 
period  which  depends  on  the  partitioning  of  the  energy  transfer  between  the  paths 
for  baroclinic  and  horizontal  shear  instabilities,  although  a  higher  frequency 
similar  to  that  for  the  baroclinic  instability  tends  to  be  superimposed  (compare 
Figs.  25b  and  11a  with  mixed  stability  to  Figs.  12c  and  13a  with  baroclinic 
instability).  In  the  mixed  stability  regime  the  primary  period  falls  between  those 
for  pure  baroclinic  and  pure  horizontal  shear  instability.  These  ideas  will  be  the 
subject  of  further  investigation  in  future  work. 

Both  the  eddy  shedding  period  and  the  eddy  diameter  for  the  reduced  gravity 
model  are  closer  to  those  for  the  bottom  topography  experiment  than  to  those  for 
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the  flat  bottom  experiment.  Since  inclusion  of  Gulf  topography  significantly 
reduces  the  transfer  of  energy  to  the  deep  water  compared  to  the  flat  bottom  model, 
our  use  of  the  reduced  gravity  model,  which  totally  prevents  such  a  transfer,  is 
both  useful  and  in  some  ways  more  realistic  than  the  use  of  a  two-layer  flat  bottom 
model.  Thus  the  reduced  gravity  model  may  be  more  representative  of  the  actual 
Gulf  than  the  two-layer  flat  bottom  model. 

Bottom  topography  in  conjunction  with  deep  water  inflow  through  the  Yucatan 
Straits  can  have  an  additional  important  effect  on  the  path  and  stability  of  the 
Loop  Current.  We  will  indicate  the  transport  distribution  (in  Sv)  between  the  two 
layers  by  the  notation  Tj/T2,  where  Ti  represents  the  upper  layer  transport  on 
inflow  and  T2  represents  the  lower  layer  transport  on  inflow.  In  the  next 
experiment  we  utilize  the  model  with  bottom  topography  described  above.  At  the  end 
of  year  6  we  gradually  increase  the  deep  water  transport  from  0  to  10  Sv  with  a 
time  constant  of  30  days.  Figure  26a  shows  the  domain-averaged  kinetic  energy  for 
years  4-8.  Note  that  the  eddy  shedding  totally  ceases  when  the  deep  transport  is 
increased  even  though  the  maximum  vertical  shear  decreases  by  only  about  10%.  The 
nearly  steady  flow  of  the  Loop  Current  is  indicated  by  the  PHA  at  day  2880  in 
Figure  26b.  When  the  deep  transport  is  again  reduced  to  zero  the  eddy  shedding 
resumes,  as  we  will  demonstrate  shortly. 

Initially,  two  hypotheses  were  considered  for  explaining  these  model  results. 
The  first  is  that  the  topographic  beta  effect  acts  to  oppose  planetary  vorticity 
advection  on  inflow  due  to  the  downward  sloping  topography.  This  would  inhibit 
westward  spreading  and  subsequently  eddy  shedding.  The  second  is  that  the  Campeche 
Bank  blocks  the  westward  spreading  of  the  Loop  and  also  inhibits  the  eddy  shedding 
process.  One  experiment  to  test  these  hypotheses  involved  totally  removing  both 
the  Campeche  Bank  and  the  bottom  slope  on  inflow.  Figure  27a  shows  the  modified 
topography  used  for  this  experiment.  Using  this  topography  we  repeated  the 
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previous  experiment  which  produced  a  steady  flow  with  a  25/10  transport  profile. 
After  three  years  of  integration  the  model  reached  a  steady  state  very  similar  to 
that  of  the  full  bottom  topography  experiment  (Fig.  26b).  Clearly,  neither  the 
bottom  slope  on  inflow  nor  the  Canpeche  Bank  topography  are  responsible  for  this 
steady  regime.  At  the  beginning  of  year  4  we  gradually  reduced  the  deep  water 
transport  to  zero  using  a  30-day  time  constant.  Eddy  shedding  resumed  as  in  the 
full  topography  case.  This  is  shown  by  the  domain  averaged  kinetic  energy  in 
Figure  27b.  In  another  experiment  we  used  the  topography  of  Figure  27a,  except 
that  we  moved  eastward  by  100  km  all  of  the  topographic  features  along  the  eastern 
boundary,  unlike  the  previous  case,  the  25/10  transport  profile  produced  eddy 
shedding.  This  demonstrates  the  importance  of  the  location  of  the  Florida  Shelf  in 
controlling  the  dynamics. 

To  elucidate  the  dynamics  of  the  transition  from  stable  to  unstable  regimes, 
we  first  examined  the  kinematics  of  the  flow  using  the  upper  layer  continuity 
equation,  (2).  In  the  steady  state  the  divergence  term,  h^V.v^,  is  balanced  by  the 
advective  term,  v^.Vh^.  In  Figure  28  we  have  plotted  the  bottom  topography,  the 
divergence  term  and  the  advective  term  for  the  steady  regime  over  a  500  x  500  km 
region  near  the  outflow  boundary.  The  large  divergence  over  the  topography 
(0.024  cm  sec-1)  is  clearly  shewn  in  Figure  28b. 

We  can  gain  further  insight  into  the  role  of  the  advective  term  in  (2)  by 
assuning  quasi-geos trophy  such  that 

V7hi  3  Vjg-Vhj  =  £  J  (h2  +  D.hj) 

■+ 

where  v^  is  the  geostrophic  velocity  and  D  is  the  height  of  the  topography  above  a 

reference  level.  The  Jacobian  may  also  be  rewritten  as  J  (n2,ni),  where  nj  is  the 

deviation  of  the  free  surface  from  its  rest  value  and  is  the  interface  deviation 
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from  its  rest  value.  Figure  28d  is  a  contour  map  of  fj  (n2,  nj) .  It  is  virtually 
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identical  to  Figure  28c,  confirming  quasi-geostrophy. 


Since  the  flow  is  quasi-geostrophic,  v^  .vh^  =  v^g  .  vhi  *  V2g.vhi.  We  will 
use  this  information  to  show  how  the  flow  in  the  lower  layer  can  prevent  the 
northward  penetration  of  the  Loop  Current.  The  magnitude  of  V2g.7hi  is  greatest 
when  relatively  strong  lower  layer  currents  flow  at  large  angles  to  contours  of 
upper  layer  thickness.  Comparison  of  Figures  26b  (contours  of  hj  -Hi)  and  26c 
(contours  of  lower  layer  pressure)  shows  that  this  occurs  when  the  Loop  Current 
approaches  the  Florida  Shelf,  since  a  current  in  the  lower  layer  tends  to  follow 
the  f/h  contours  of  the  shelf.  If  this  current  is  strong  enough  for  the  advection 
to  balance  the  divergence  associated  with  the  approaching  Loop  in  (2) ,  then  the 
interface  deepening  and  Loop  Current  penetration  in  this  region  are  halted. 

How  does  this  affect  the  penetration  of  the  remainder  of  the  Loop?  We  have 
seen  two  cases  (Figs.  26  and  27)  where  the  Loop  penetration  ceased  entirely.  In  a 
third  case,  where  the  Florida  Shelf  had  been  moved  further  eastward,  the  Loop 
Current  penetrated  far  into  the  Gulf  and  shed  eddies.  The  Florida  Shelf  in 
conjunction  with  the  lower  layer  flow  effectively  changes  the  port  separation  by 
locally  limiting  the  northward  penetration  of  the  current.  If  this  results  in 
Rq  >  2  when  we  measure  2  Lp  as  the  distance  between  the  center  of  the  inflow  port 
and  the  point  of  maximum  v^.Vhi  (Fig.  28) ,  then  the  current  falls  into  regime  N 
with  source-sink  flow,  a  regime  discussed  in  section  4a  (Fig.  18) .  This  is 
illustrated  by  comparing  Figs.  17  and  26b.  Figure  17  is  a  reduced  gravity  case 
with  a  land  mass  inserted  which  was  designed  to  match  the  location  of  the  Florida 
Shelf. 

At  the  beginning  of  year  4  we  reduced  the  deep  water  inflow  transport  in  the 
experiment  shown  in  Figure  27.  Soon  after,  the  Loop  Current  penetrated  northward 
into  the  Gulf,  spread  westward,  and  proceeded  to  shed  eddies  in  a  regular  fashion. 
When  the  lower  layer  inflow  transport  was  reduced,  the  current  in  the  lower  layer 
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along  the  Florida  Shelf  weakened  and  v^. vh^  was  no  longer  able  to  balance  the 
divergence  term.  Thus,  the  interface  deepened  and  the  Loop  Current  penetrated 
northward. 

Figure  26a  shows  a  case  where  the  Loop  Current  had  already  penetrated  far  into 
the  Gulf  before  the  lower  layer  inflow  through  the  Yucatan  Straits  was  increased. 
When  this  increase  occurred,  a  current  in  the  lower  layer  developed  along  the 
Florida  Shelf  slope  following  f/h  contours  and  the  advection  (vi- vh^)  began  to 
exceed  the  divergence  term  in  the  continuity  equation  where  the  Loop  Current 
crossed  the  shelf  slope,  causing  the  interface  to  shallow  there.  Thus,  the  Loop 
Current  began  to  retreat  southward  until  an  equilibrium  occurred  near  the  southern 
end  of  the  shelf. 

It  is  also  of  interest  to  investigate  the  changes  in  the  vorticity  balance  due 
to  topography.  We  find  same  notable  differences  from  the  results  of  Holland 
(1973) .  From  a  vorticity  analysis  we  can  obtain  term  balances  for  each  layer  and 
for  the  total  fluid  coluim.  The  steady  state  mass  transport  vorticity  equations 
for  each  layer  and  for  the  total  depth  are 

(NL)j  +  BVj  +  g  J(npn2)  -  A  v2^  =  0  (21) 

(NL)2  +  ev2  +  (g-g1)  JUj.D)  -  (g-g1)  J( n x ,n2)  +  g'J(n2,D)  -  A  v2^  =  0  (22) 

p i (NL. ) i  +  p 2 ( NL ) 2  +  8(p1V1  +  p2V2)  +  p2(g-g')  0(n1,D)  +  p2g'J(n2,D)  (23) 

-  A  v2  (p^1  +  p2c2)  =  0 

We  have  represented  the  sum  of  the  nonlinear  terms  in  the  vorticity  balance  by  the 
letters  NL,  and  here  =  -  Uiy.  Note  that  in  each  equation  there  are 

terms  representing  pressure  torques.  For  the  total  mass  transport,  (23) ,  the 
bottom  pressure  torque  is  represented  by  two  terms.  The  first  term  is  analogous  to 
Holland's  (1973)  surface  pressure  torque  term  and  the  second  corresponds  to  his 
density  compensation  term. 
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Figure  29  provides  a  plot  of  the  term  balances  in  Eq.  21-23  along  a  section 
crossing  the  Loop  Current  from  x  *  600  km  to  x  *  1600  km  at  y  »  197.5  km.  This 
section  passes  directly  over  the  area  of  strong  divergence  shown  in  Figure  28. 

For  both  the  deep  layer  and  for  the  sum  of  the  layers  the  individual  pressure 
torque  terms  have  been  combined  for  plotting  purposes. 

In  the  upper  layer,  across  that  part  of  the  stream  flowing  over  steep 
topography,  the  "thermocline  pressure  torque"  very  nearly  balances  the  nonlinear 
terms.  Planetary  vorticity  advection  is  surprisingly  unimportant,  as  is  vorticity 
diffusion.  In  the  deep  water  the  length  scales  are  considerably  shorter  than  at 
the  surface,  being  dominated  by  the  scale  of  the  topographic  slope.  The  term 
balance  is  not  as  clear-cut  in  the  deep  water,  although  the  pressure  torques  still 
tend  to  balance  the  nonlinear  terms.  The  sum  of  mass  transport  vorticity  from  each 
layer  yields  a  clear  balance  between  the  bottom  pressure  torque  and  the  nonlinear 
terms  over  the  topography.  This  is  in  marked  contrast  to  Holland  (1973)  who  found 
that,  in  the  Gulf  Stream  over  the  continental  slope,  the  bottom  pressure  torque  was 
balanced  by  planetary  vorticity  advection.  The  nonlinear  terms  in  his  experiment 
contributed  less  than  one  percent  to  the  total. 

In  terms  of  vorticity,  the  dynamics  for  this  regime  of  the  Loop  Current,  when 
deep  transpert  is  significant,  may  be  sunmariEed  as  follows.  The  Loop  Current 
becomes  "locked"  into  the  steady  regime  when  the  deep-water  transport  is 
sufficiently  large  that  it  produces  strong  pressure  torques  which  balance  relative 
vorticity  advection  over  the  steep  slope,  in  that  case  deepening  of  the  layer 
interface  is  prevented  and  the  Loop  Current  is  unable  to  penetrate  into  the  Gulf  - 
a  prerequisite  for  spreading  westward  into  the  unstable  configuration  which  leads 
to  eddy  shedding.  When  the  deep  transport  is  reduced  and  the  pressure  torques 
decrease,  the  steady  balance  breaks  down  and  the  interface  begins  to  deepon.  Once 
this  deepening  process  begins,  the  Loop  is  able  to  penetrate  northward  into  the 
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Gulf,  spread  westward  under  the  influence  of  Rossby  wave  dynamics,  and  reach  an 
unstable  configuration. 

Our  results  strongly  suggest  that  the  deep  water  transport  through  the  Yucatan 
Straits  may  be  of  crucial  ijqportance  for  the  control  of  Loop  Current  penetration 
and  stability.  Rather  than  monitoring  the  transport  and  velocities  in  the  upper 
part  of  the  water  colunn,  it  may  be  more  important  to  monitor  the  deep  transport. 
Significant  changes  in  deep  water  transport  through  the  Yucatan  Straits  may  be 
reflected  by  marked  changes  in  the  behavior  of  the  Loop  Current.  The  recently 
reported  observations  of  southward  flow  near  the  bottom  of  the  Yucatan  Channel 
during  a  one-month  period  (Hansen  and  Molinari,  1979)  indicate  the  need  for 
additional  modeling  and  field  studies  of  the  deep  flow  and  its  influence  on  the 
Loop  Current. 

c.  Time-dependent  inflow 

Our  investigations  have  demonstrated  that  the  growth  and  instability  of  the 
Loop  Current  do  not  require  time-dependent  inflows  as  an  essential  element  of  the 
dynamics.  Nevertheless,  it  is  possible  that  the  variability  of  total  transport 
and/or  current  profiles  through  the  Yucatan  Straits  significantly  influences  the 
actual  Loop  Current  behavior.  Unfortunately,  there  exist  few  direct  current 
measurements  in  the  Yucatan  Straits  to  provide  estimates  of  spatial  or  tenporal 
variability  of  the  inflow  on  time  scales  comparable  to  the  eddy  shedding  period. 
Existing  geos trophic  transport  calculations  based  on  hydrographic  sections  suffer 
from  the  classical  problem  of  choice  of  reference  level.  Many  investigators, 
including  Schlitz  (1973),  Morrison  and  Nowlin  (1977),  Molinari  and  Yager  (  1977), 
Maul  (1977,  1979),  and  Hansen  and  Molinari  (1979)  have  estimated  geostrophic 
transport  through  the  Yucatan  Straits  on  time  scales  from  days  to  months.  Molinari 
and  Yager  (1977)  found  from  repeated  hydrographic  surveys  across  the  Yucatan 
Straits  during  a  26-day  period  during  May  1972,  that  the  mean  transport  relative  to 


700  m  was  18.2  (+  1.3)  Sv  with  a  7%  variation  in  transport  over  the  period.  Hansen 
and  Molinari  (1979) ,  using  both  hydrographic  data  and  direct  current  meter  data, 
estimated  total  transport  across  the  Straits  to  be  25.8  Sv  during  October  3-5,  1970 
and  27.1  Sv  during  October  26-27,  1970.  The  5%  difference  in  transport  was 
considered  to  be  within  observational  error.  The  cross-stream  position  of  the 
maximum  current  above  700  m  varied  from  near  the  center  of  the  stream  to  near  the 
western  boundary  during  the  three-week  period.  The  1972-73  estimates  of 
geostrophic  transport  relative  to  700  m  by  Maul  (1979)  indicated  highest  transports 
and  current  speeds  in  December,  February,  and  March  and  minimum  values  in  May,  June 
and  July.  His  finding  is  counter  to  the  cycle  reported  by  Fuglister  (1951)  and 
Cochrane  (1965)  who  reported  speed  maxima  during  the  surrmer  months. 

Although  there  are  no  statistically-signif icant,  annual  time-scale 
measurements  of  transport  through  the  Yucatan  Straits,  proxy  data  from  the  Florida 
Straits  obtained  by  Niiler  and  Richardson  (1973)  do  exist.  Seven  years  of  directly 
measured  transport  and  horizontal  current  data  from  Miami  to  Bimini  were  fitted  to 
an  annual  sinusoid,  giving  a  mean  total  transport  of  29.5  Sv  with  an  amplitude  of 
4.1  Sv  and  a  peak  transport  in  early  June.  Importantly,  there  was  as  much  variability 
within  a  particular  season  as  from  season  to  season. 

For  the  Yucatan  Straits,  the  oceanographic  observations  are  presently 
inadequate  to  provide  unbiased  estimates  of  cross-stream  current  profiles  on  the 
eddy  shedding  time  scale.  At  best  the  data  give  short  term  estimates  of  mean  total 
transport,  cross-stream  transport  profiles  and  transport  fluctuations.  The  data 
suggests  that  both  cross-stream  velocity  profiles  and  total  transport  have  short 
term  (days  to  weeks)  fluctuations  comparable  to  the  seasonal  and  longer-term 
fluctuations. 

To  test  the  hypothesis  that  Loop  Current  growth  and  eddy  shedding  are 
correlated  with  transport  fluctuations  through  the  Yucatan  Straits  on  the  seasonal 
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time  scale,  we  have  conducted  numerical  experiments  in  which  a  time-varying 
component  was  super  imposed  on  the  mean  inflow  transport.  Unlike  Wert  and  Reid 
(1972) ,  who  varied  the  velocity  profile  markedly  over  an  annual  cycle  while  leaving 
the  total  transport  constant,  we  have  not  considered  in  detail  time-dependent 
fluctuations  in  the  velocity  profile  across  the  Straits.  However,  as  discussed  in 
section  3,  it  is  the  speed  of  the  Loop  Current  rather  than  its  relative  vorticity 
distribution  which  is  most  important  in  determining  the  eddy  shedding  rate,  the 
diameter  of  the  eddies,  and  the  eddy  intensity.  For  example,  the  standard  20  Sv 
reduced  gravity  experiment  with  a  parabolic  inflow  profile  and  maximum  inflow  speed 
of  72.9  cm  sec~l  yielded  an  eddy  shedding  period  of  327  days  and  a  characteristic 
eddy  diameter  of  373  km.  The  identical  20  Sv  experiment  with  a  uniform 

cross-stream  profile  and  an  inflow  speed  of  63  an  sec-1  was  also  integrated  to  statistical 
equilibrium.  The  conservation  of  potential  vorticity  dynamics  discussed  in  section 
3  indicate  that  the  eddy  diameter  and  eddy  shedding  period  for  this  flat  profile 
case  should  be  smaller  by  the  square  root  of  the  ratio  of  the  maximum  inflow  speeds 
for  the  two  cases.  This  yields  a  shedding  period  for  the  flat  profile  case  of  304 
days  and  an  eddy  diameter  of  347  km.  The  observed  period  for  the  flat  profile 
experiment  is  308  days,  the  observed  diameter  348  km.  Thus  we  are  able  to  explain 
the  results  from  these  two  experiments  with  differing  cross-stream  profiles  on 
inflow  with  only  a  knowledge  of  the  maximum  current  speeds  in  each  profile.  These 
results  suggest  that  time  variations  in  inflow  speed  at  the  core  of  the  current  are 
considerably  more  important  to  the  shedding  process  than  time  variations  of  the 
detailed  shape  of  the  cross-stream  inflow  profile.  We  shall  concentrate  on  this 
aspect  of  the  time-dependent  forcing.  The  experiments  described  in  this  subsection 
use  a  time-varying  total  inflow  transport  with  a  parabolic  profile  for  V  and  relatively 
small  time-dependent  variations  in  the  shape  of  the  cross-stream  velocity  profile 
associated  with  upwelling  and  downwelling  at  the  inflow  port. 

Four  experiments  were  performed  in  which  a  sinusoidally-varying  component  was 
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superimposed  on  the  mean  20  Sv  inflow  of  the  standard  reduced  gravity  model. 
Peak-to-peak  amplitude  of  the  time-dependent  component  was  10  Sv.  This  is  at  least 
as  large  as  the  observed  maximum  transport  fluctuation  observed  through  the  Yucatan 
Straits  on  any  time  scale.  Each  of  the  four  experiments  was  run  to  statistical 
equilibrium  with  a  period  of  360,  240,  180,  or  110  days  for  the  time  varying 
component  of  the  inflow.  Figure  30  summarizes  the  results  of  these  experiments  in 
terms  of  potential  energy  density  plots  for  the  standard  case  with  constant  inflow 
transport  and  for  the  360,  240  and  180  day  period  time-dependent  forcing.  These 
curves  were  calculated  from  time  series  of  PHA  constructed  at  a  point  in  the  model 
where  shed  eddies  were  always  observed  to  pass.  In  these  experiments  the  point  was 
located  at  x  =  810  km,  y  =  421.9  km. 

The  327  day  eddy  shedding  period  for  the  constant  inflow  case  shows  up  clearly  as 
a  pronounced  maximum  in  the  spectrum  (Fig.  30a).  Harmonics  near  periods  of  164,  109, 
and  82  days  are  also  clearly  evident.  We  will  term  the  327  day  period  the  "natural" 
period  of  oscillation  for  the  standard  case.  When  we  add  the  annual  forcing  to  the 
mean  inflow,  the  maximum  spectral  peak  shifts  to  the  one-year  period,  as  shown  in 
Figure  30b.  There  is  no  peak  at  the  natural  period.  Harmonics  at  6,  4  and  3 
months  are  clearly  evident.  The  annual  cycle  for  the  eddy  shedding  is  also 
observed  when  individual  eddies  are  monitored  in  their  separation  from  the  Loop 
Current.  This  suggests  that  the  maximum  in  the  spectrum  may  consistently  appear  at 
the  forcing  frequency.  However,  in  the  eight  month  forcing  experiment  no  dominant 
spectral  maximum  at  the  forcing  frequency  is  found  (Fig.  30c) .  Rather,  we  observe  a 
broad  maximun  centered  near  the  natural  frequency.  The  eddy  shedding  periods  we 
actually  observe  in  this  model  realization  range  fran  7  to  13  months,  with  considerable 
variability  from  one  eddy  cycle  to  the  next.  This  is  shown  in  the  phase  plots  for 
the  PHA  in  Figure  31a.  This  range  of  variability  for  the  eddy  shedding  is 
comparable  to  that  observed  (Behringer,  et  al.  1977). 
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The  experiment  with  semi-annual  forcing  produced  a  striking  result.  The 
eddy  shedding  period  was  nearly  annual  (Fig.  30d) .  As  discussed  in  section  3,  this  is 
because  the  time  scale  for  the  growth  of  the  horizontal  shear  instability  is  short 
compared  to  the  time  required  for  the  loop  Current  to  penetrate  into  the  Gulf  and  bend 
westward  into  an  unstable  configuration.  In  this  experiment  the  first  maximum  in 
inflow  transport  for  the  year  occurred  when  the  Loop  was  just  beginning  to 
penetrate  into  the  Gulf.  The  transport  maximum  and  subsequent  decrease  did  not 
coincide  with  an  unstable  configuration.  However,  by  the  time  the  second  yearly 
maximum  in  inflow  transport  occurred,  the  Loop  had  penetrated  well  into  the  Gulf 
and  the  time  variation  in  the  inflow  was  able  to  exert  some  influence  on  the  timing 
of  the  eddy  shedding. 

Because  the  observed  period  is  356  days  while  the  forcing  period  is  180  days, 
the  forcing  and  the  shedding  gradually  shift  out  of  phase  until  the  first  transport 
maximum  of  the  year  rather  than  the  second  influences  the  timing  of  the  eddy 
shedding.  Figure  31b  shows  the  phase  plot  of  PHA  for  years  5  through  15  for  this 
experiment.  Note  that  during  years  5  through  10  the  eddy  passes  through  our 
north-south  section  with  maximum  amplitude  slightly  after  the  beginning  of  a  year. 

The  time  between  eddy  sheddings  is  356  days.  This  yields  a  very  sharp  peak  in  the 
potential  energy  spectrum,  as  shown  in  Figure  30d.  However,  by  year  11  we  see  a 
transition  in  the  time  of  eddy  passage  in  Figure  31b,  and  by  year  12  the  eddies  are 
passing  through  our  section  with  maximum  amplitude  just  after  the  middle  of  a  year. 
Following  the  transition  the  eddy  shedding  period  is  again  approximately  1  year. 

Were  we  to  recalculate  our  spectrum  for  this  experiment  using  a  very  long  record 
with  equal  numbers  of  shed  eddies  for  each  phase  of  eddy  shedding,  the  maximum  at 
the  annual  period  would  be  reduced  and  a  semi-annual  peak  should  become  apparent. 

The  final  experiment  with  a  110  day  forcing  period  yielded  results  similar  to 
the  semi-annual  forcing  except  that  the  strong  spectral  rraxinam  occurred  at  the  330  day 


period  with  higher  harmonics  also  containing  significant  energy. 

These  experiments  indicate  that  time-dependent  fluctuations  with  periods  very 
near  the  natural  eddy  shedding  period  may  shift  the  eddy  shedding  to  nearly  the 
forcing  period.  For  inflow  forcing  with  periods  significantly  different  from  the 
natural  period  or  its  harmonics,  the  eddy  shedding  may  become  erratic,  with 
considerable  variation  in  the  shedding  period.  The  short  time  scale  for  growth  of 
the  Loop  Current  instability  compared  to  the  time  scale  for  penetration  of  the  Loop 
into  the  Gulf  has  significant  implications  for  Loop  Current  response  to  seasonal  time 
scale  fluctuations  in  inflow  forcing.  If  the  loop  is  close  to  an  unstable  configuration, 
transport  fluctuation  at  inflow  may  substantially  influence  the  timing  for  eddy 
break-off. 

It  should  be  noted  that  our  time-dependent  experiments  were  conducted  with  the 
reduced  gravity  model.  However,  from  section  4b  it  is  clear  that  changes  in  deep 
water  transport  through  the  Yucatan  Straits  may  lead  to  cessation  of  Loop  Current 
penetration  and  eddy  shedding.  Thus,  relatively  small  fluctuations  in  deep  water 
transport  on  seasonal  or  annual  time  scales  could  modulate  the  Loop  Current 
penetration  and  eddy  shedding. 
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5.  Summary  and  conclusions 


We  have  used  three  nonlinear  numerical  models  to  study  the  dynamics  of  the 
Loop  Current-eddy  system  in  the  Gulf  of  Mexico:  two-layer,  barotropic  and  reduced 
gravity.  A  second  active  layer  allowed  the  possibility  of  baroclinic  instability, 
and  allowed  inclusion  of  both  the  pycnocline  and  idealized  topography  of  the  Gulf. 
The  barotropic  and  reduced  gravity  models  were  used  to  demonstrate  the  individual 
behavior  of  the  external  and  internal  modes,  and  to  provide  insight  on  how  they 
interact  in  the  two-layer  model.  They  also  allowed  us  to  study  certain  dynamical 
features  of  the  system  in  the  simplest  context.  Because  of  the  economy  of  the 
semi-implicit  free  surface  models,  it  was  possible  to  perform  over  100  multi-year 
experiments  to  investigate  the  dynamics  of  the  Loop  Current-eddy  system.  Typically 
the  models  were  integrated  3  to  5  years  to  statistical  equilibrum  on  a 
1600  x  900  km  domain  with  resolution  20  x  18.75  km.  Prescribed  inflow  through  the 
Yucatan  Channel  was  compensated  by  outflow  through  the  Florida  Straits. 

A  longstanding  hypothesis  is  that  the  Loop  Current  sheds  eddies  in  response  to 
quasi-annual  variations  in  the  inflow  through  the  Yucatan  Straits.  However,  we 
first  searched  for  purely  internal  mechanisms  to  produce  the  eddy  shedding. 

Instead  of  starting  with  the  simplest  possible  model  and  progressing  toward  more 
complicated  experiments,  we  designed  our  first  experiment  to  be  as  realistic  as 
possible  within  the  context  of  the  two-layer  model.  This  meant  including  idealized 
topography  of  the  Gulf.  However,  in  our  search  for  purely  internal  mechanisms  for 
the  eddy  shedding,  we  used  a  time  invariant  inflow  profile.  We  found  that  the 
model  can  realistically  simulate  observed  features  of  the  eddy  shedding  by  the  Loop 
Current  when  model  parameters  realistic  for  the  Gulf  of  Mexico  are  used,  including 
the  size,  amplitude,  and  movement  of  the  eddies.  We  also  found  that,  contrary  to 
an  earlier  hypothesis,  time  variations  in  the  inflow  are  not  required  for  the  eddy 
shedding  to  occur  at  a  realistic  rate. 


63 


We  next  showed  that  the  reduced  gravity  model  is  the  simplest  model  which  can 


simulate  the  basic  features  of  the  eddy  shedding  cycle  of  the  Loop  Current.  Thus, 
baroclinic  instability  is  not  an  essential  element  of  the  dynamics.  Instead,  there 
is  a  horizontal  shear  instability  of  the  first  internal  mode.  We  also  found  that 
differential  rotation  is  essential  for  eddy  shedding  to  occur  in  a  realistic 
manner.  When  we  set  £=0,  the  solution  evolved  to  a  steady  source-sink  flow  with  no 
eddy  shedding  or  westward  spreading.  This  was  true  even  when  we  augmented  the 
inflow  with  a  large  westward  component  and  an  annual  variation  larger  than 
observed. 

Fortunately,  we  were  able  to  perform  enough  reduced  gravity  experiments  to 
test  a  number  of  dynamical  hypotheses  using  regression  analysis.  We  found  that  the 
westward  speed  of  the  shed  eddies  was  correlated  with  the  nondispersive  internal 
Rossby  wave  speed  with  a  correlation  coefficient  of  0.991,  but  it  was  necessary  to 
include  dispersion  to  obtain  quantitative  agreement  of  the  speeds.  Theoretical 
expressions  for  the  eddy  diameter  and  the  penetration  distance  of  the  Loop  Current 
into  the  Gulf  were  obtained  by  extending  a  potential  vorticity  argument  on  a 
3-olane  by  Reid  (1972).  The  theoretical  diameter,  2r,  was  consistent  with  beta 
Rossby  number,  Rg  =  Vjn/tBr2)  =  1,  and  predicted  the  diameter  of  the  eddies  found 
numerically  with  an  average  error  of  6%,  based  only  on  a  knowledge  of  the  value  of  0 
and  the  maximum  speed  of  the  current  at  inflow,  Vin.  A  simplification  of  the 
expression  for  penetration  distance  was  tested  using  regression.  We  obtained  an 
average  fit  error  of  2%  and  a  linear  correlation  of  0.986,  using  eddy  diameters 
obtained  numerically.  Hence,  we  have  demonstrated  the  importance  of  differential 
rotation  in  determining  the  diameter  of  the  eddies,  and  the  penetration  of  the 
model  Loop  Current  into  the  Gulf. 

Two  time  scales  are  associated  with  the  eddy  shedding:  the  time  for  the  Loop 
Current  to  penetrate  into  the  Gulf  and  spread  westward  into  an  unstable 
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configuration,  and  the  much  shorter  time  scale  for  the  growth  of  a  horizontal  shear 
instability  as  the  eddy  separates  from  the  Loop  Current.  The  Loop  Current  does  not 
necessarily  meet  any  criterion  for  instability  over  much  of  an  eddy  cycle.  The 
eddy  shedding  period  is  a  small  multiple  of  the  time  required  for  an  eddy  to  move 
the  diameter  of  the  eddy  after  shedding.  This  multiple  depends  on  the  inflow  angle 
of  the  current.  Thus,  in  the  most  realistic  regime,  the  eddy  shedding  rate  depends 
on  the  internal  Rossby  wave  speed,  an  eddy  diameter  derived  from  conservation  of 
potential  vorticity  on  a  8-plane,  and  the  angle  of  the  inflow.  To  a  lesser  extent 
it  also  depends  on  the  Reynolds  number.  As  the  flow  becomes  more  inertial,  the 
eddy  shedding  period  is  increased,  perhaps  a  counter-intuitive  result.  This  is  due 
to  the  way  friction  modifies  the  results  predicted  by  conservation  of  potential 
vorticity  on  a  6-plane.  Correlations  of  0.95  to  0.98  and  mean  fit  errors  of  6  to 
8%  were  obtained  in  our  regression  predictions  of  eddy  shedding  period. 

The  eddy  shedding  can  be  prevented  by  reducing  the  Reynolds  number 
sufficiently.  However,  the  Loop  Current  still  spreads  far  westward.  The  steady 
state  solution  for  a  highly  viscous  case  was  found  to  be  almost  the  same  as  the 
mean  over  an  eddy  cycle  of  a  lower  viscosity  case  which  shed  discrete  eddies  of 
large  amplitude.  In  a  few  of  the  two-layer  experiments,  baroclinic  and  mixed 
instabilities  were  encountered,  but  experiments  dominated  by  a  horizontal  shear 
instability  of  the  internal  mode  produced  the  most  realistic  results.  For 
sufficiently  high  Reynolds  numbers  the  shear  instability  occurred  in  both  the 
barotropic  and  reduced  gravity  models.  However,  for  realistic  parameter  values 
eddy  shedding  occurred  in  the  two-layer  and  reduced  gravity  models,  but  not  in  the 
barotropic  model. 

Experiments  using  the  reduced  gravity  model  showed  that  while  the  solution 
depends  on  the  maximum  velocity  at  inflow,  it  is  relatively  insensitive  to  the 
shape  of  the  inflow  profile  or  to  the  location  of  the  ports.  When  the  ports  were 
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relocated  by  several  hundred  kilometers,  the  Loop  penetration,  the  eddy  shedding 
period,  and  the  eddy  diameter,  amplitude  and  movement  were  little  affected,  so  long 
as  the  western  boundary  did  not  interfere  with  the  shedding  process,  and  the  ports 
were  not  unrealistically  close.  The  solution  was  also  quite  insensitive  to  the 
width  of  the  outflow  port.  Even  when  the  entire  eastern  boundary  was  opened,  the 
outflow  was  confined  to  a  current  adjacent  to  the  southern  boundary,  suggesting  a 
dynamical  rather  than  a  geological  explanation  for  the  location  of  the  Florida 

Straits.  These  results  are  consistent  with  the  domination  of  the  model  behavior  by 
conservation  of  potential  vorticity  dynamics. 

Scxne  of  the  diverse  results  from  the  barotropic  flat  bottom  model  and  the 
reduced  gravity  model  have  been  summarized  in  a  regime  diagram  with  Reynolds 
number,  Re  vs.  beta  Rossby  number,  Rg  (Fig.  18) .  The  length  scale  for  Re  is  one 
half  the  port  width,  for  %  one  half  the  port  separation.  For  low  Re  and  % 

(regime  w)  the  Loop  Current  spreads  westward  and  evolves  to  a  steady  state  without 
shedding  eddies.  For  high  Re  and  low  Rg  (regime  E)  the  Loop  Current  sheds  eddies 
with  no  time  variation  in  the  external  forcing.  Near  the  transition  between 
regimes  E  and  W  the  eddies  are  superimposed  on  a  westward  spreading  loop.  For  high 
Rg  (regime  N)  the  Loop  Current  evolves  to  a  steady  source-sink  flow  with  no 
westward  spreading  or  eddy  shedding.  The  f-plane  solution  belongs  to  this  regime. 
The  port  separation,  2Lp,  is  approximately  Lp<r/'v/2  for  this  regime  to  occur,  where 

j, 

r  =  (Vc/8)  2 ,  a  theoretical  eddy  radius,  and  vc  is  the  speed  at  the  core  of  the 
current. 

The  Florida  Shelf  topography  plays  a  crucial  role  in  the  Loop  Current-eddy 
shedding  dynamics  in  certain  realistic  regions  of  parameter  space.  If  the  lower 
layer  inflow  through  the  Yucatan  Straits  is  sufficiently  high,  then  the  Loop 
Current  evolves  to  a  steady  source-sink  flow  resembling  the  N  regime  with  no 
westward  spreading  or  eddy  shedding.  The  current  in  the  lower  layer  tends  to 
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follow  the  f/h  contours  of  the  Florida  Shelf  and  intersects  the  Loop  Current  in  the 
upper  layer  at  large  angles.  A  kinematic  analysis  shows  that  in  this  situation  the 
advection  term  can  balance  the  divergence  term  in  the  continuity  equation  for  the 
layered  model,  if  the  current  in  the  lower  layer  is  sufficiently  strong.  Locally, 
this  prevents  the  deepening  of  the  interface  associated  with  the  northward 
penetration  of  the  Loop  Current.  This  effectively  reduces  the  port  separation  and 
forces  the  current  into  the  N  regime.  When  the  Florida  Shelf  topography  is  moved 
100  km  eastward,  the  westward  spreading  and  the  eddy  shedding  reappears.  In  terms 
of  vorticity  dynamics,  the  northward  penetration  of  the  Loop  Current  is  halted  when 
the  interaction  between  the  topography  of  the  Florida  Shelf  and  the  pressure  field 
results  in  a  near  balance  between  the  pressure  torques  and  the  nonlinear  terms  in 
the  mass  transport  vorticity  equation. 

Another  important  role  of  the  topography  is  to  inhibit  baroclinic  instability. 
As  a  result  numerical  experiments  using  the  two- layer  model  with  idealized  Gulf 
topography  resemble  those  using  the  reduced  gravity  model  more  closely  than  those 
using  the  two-layer  flat  bottom  model,  although  the  movement  of  the  shed  eddies  is 
significantly  modified  by  the  introduction  of  topography. 

The  response  of  the  current  to  time-varying  inflows  was  ccnplex.  Importantly, 
the  eddy  shedding  frequency  was  not  dominated  by  the  forcing  frequency  for 
realistic  transport  fluctuations.  The  natural  frquency  was  of  much  greater 
importance  in  determining  the  eddy  shedding  rate  in  these  cases.  For  example,  the 
reduced  gravity  model  with  a  constant  20  Sverdrup  inflow  yielded  a  natural  eddy 
shedding  period  of  327  days.  When  a  semi-annual  sine-wave  fluctuation  with  a  10  Sv 
peak-bo-peak  amplitude  was  super  imposed  on  the  mean  inflow,  the  eddy  shedding 
period  was  356  days. 
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APPENDIX  A 


List  of  Symbols  for  Model  Equations 


A 

D 

f 

g 

g’ 

hl»h2 

Hi,H2(x,y) 

t 

“l#u2 

vl»v2 


horizontal  eddy  viscosity 

height  of  bottom  topography  above  a  reference  level 

Coriolis  parameter 

acceleration  due  to  gravity 

reduced  gravity,  g  [p2  -  Pi)/p 

instantaneous  local  thickness  of  the  layers 

initial  thickness  of  the  layers 

time 

x-directed  components  of  current  velocity 
y-directed  components  of  current  velocity 


VX,V2 


x,  y,  z 


6 

At 

AX,  Ay 
nl 

n2 

P,  Pi,  p2 

iiX,  Hy 


hivi#h2v2 

tangent  plane  Cartesian  coordinates:  x  positive  eastward, 
y  positive  northward,  z  positive  upward 
differential  rotation,  df/dy 
time  increment  in  the  numerical  integration 
horizontal  grid  increments 

free  surface  anomaly;  height  of  the  free  surface  above  its 
initial  uniform  elevation;  m  hi  +  h2  -  Hi  -  H2 
n2  *  Hi  +  ni  -  hi  *  h2  -  H2  -  -  PHA 
densities  of  sea  water 

x  and  y  directed  tangential  stresses  at  the  top  (i)  and 
bottom  (i+1)  of  layer  i. 


Additional  symbols  are  defined  in  Tables  1  and  2 


APPENDIX  B 


Numerical  Design  of  the  Semi-Inplicit  Gulf  of  Mexico  Models 


The  model  equations,  (1)  through  (3) ,  are  solved  numerically  on  a  staggered 
grid.  The  smallest  representation  of  the  ”C”  stencil  (e.g. ,  Mesinger  and  Arakawa 
1976)  is 

v 

u  h  u 


Defining  the  operators, 

Sm2(VI(z))=  Jw(z+mA/2)  -  W(z-mA/2)J  , 
W(z)  =  j  |w(z+mA/2)  +  W{z-mA/2)J  , 


and 


=  6  (  )  6  ( 
xx  +  yy 


where  W  is  a  function  of  the  discrete  variable,  z,  a  refers  to  a  space  or  time 
increment,  and  m  is  an  integer,  the  finite  difference  form  of  the  model 

equations  is 

- x 

rx  -x.  -x  -v 


S2t  Ui  =  ‘5x  (Ui  ui>  -  4y  <Vi  ui>  + 


—X 


t-At 


(Al) 


--y  -x 


-y  -y 


-X 


2t 


V1  -  “6x(Ui  vi>  -  6y  <Vi  V  -  fuihi 


~2t 


-M-  Pi 


-  <ivM  V,  +  *  *6‘ 


t-At 


(A2) 


-x  -y 

where  i*l,2,  u»U/h  ,  v»V/h  ,  the  Hi  are  constants,  and  all  differences  are 
centered  in  space  and  time  unless  otherwise  noted. 


The  time  differencing  is  leap-frog,  but  an  Euler  difference  is  used  for  start 
and  restart  to  prevent  development  of  the  time  splitting  associated  with  this 
scheme.  Every  ten  days  (160  time  steps)  the  solution  for  two  adjacent  time  steps 
is  averaged,  the  solution  written  to  a  history  tape,  and  an  Euler  restart 
performed.  This  procedure  has  virtually  no  effect  on  our  numerical  solutions.  For 
advection  a  variant  of  the  Lilly  (1965)  scheme,  also  known  as  Scheme  C  from 
Grammeltvedt  (1969) ,  is  used.  The  Holland  and  Lin  (1975)  scheme  is  used  for  the 
Coriolis  force.  Horizontal  friction  is  lagged  in  time  for  computational  stability. 
The  numerical  scheme  conserves  mass  and  total  energy  except  for  the  time  average 
and  restart  procedure  and  explicit  sources  and  sinks  due  to  horizontal  friction  and 
to  inflow  and  outflow  through  the  ports. 


The  semi-implicit  treatment  of  the  external  and  internal  gravity  waves  is 
patterned  after  Kwizak  and  Robert  (1971)  and  O'Brien  and  Hurlburt  (1972)  and  allows 
us  to  greatly  exceed  the  usual  time  step  limitation  imposed  by  gravity  waves  in 
explicit  models.  By  taking  the  divergence  of  the  momentum  equations  and 
substituting  into  the  continuity  equations,  the  model  equations  can  be  solved  in 
favor  of  ni  and  h^  to  form  Helmholtz  equations  which  are  coupled  in  the  vertical: 


V  (At 


3Uj 

3X~ 


3  V 
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ay 
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,  (A5) 
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where  we  have  already  performed  the  time  differencing,  and  ai»l/gH£At2.  The  left 
hand  side  of  (A4)  and  (A5)  are  unknown  values  at  time  level  n+1,  and  the  right  hand 
side  is  the  known  forcing  function  at  time  levels  n  and  n-1.  Uj*,  Vi*,  and  Hi*  are 
the  known  forcing  functions  from  the  x -momentum,  y-momentum,  and  continuity 
equations.  To  facilitate  the  solution  of  Helmholtz  equations  (A4)  and  (A5)  in  the 
two-layer  model,  they  are  decoupled  in  the  vertical.  Use  of  Lagrange  multipliers 
is  a  convenient  way  to  do  this,  although  care  is  required  to  avoid  large  roundoff 
error  in  the  calculation  of  the  new  coefficients  of  the  Helmholtz  terms  in  the 
decoupled  equations. 


From  (A4)  and  (A5)  it  appears  that  we  must  know  the  boundary  values  for  ni  and 
or  their  normal  derivatives  at  time  levels  n-1  and  n+1  before  we  can  obtain  an 
interior  solution.  These  are  not  generally  known.  In  the  limit  as  we  approach  a 
solid  no-slip  boundary  the  normal  momentum  equation  reduces  to  a  balance  between 
the  pressure  gradient  force,  the  wind  stress,  and  the  Laplacian  horizontal 
friction.  The  situation  is  even  worse  at  the  eastern  port  where  the  normal 
velocity  at  the  boundary  of  the  model  domain  is  self-determined.  This  difficulty 
is  easily  overcome  without  further  approximation,  if  we  first  difference  the 
primitive  equations  in  space  and  time,  and  then  form  the  finite  difference  analogue 
of  the  Helmholtz  equations  algebraically.  This  procedure  also  eliminates  ambiguity 
in  defining  Ax,  Ay  on  a  grid  with  stretched  coordinates .  It  will  be  illustrated  for 
the  barotropic  model.  The  primitive  equations  can  be  written  in  finite  difference 
form  as 


1+4J  +  (h1+U  '  bl.J*  "  U*+4,j 


AX .  . 

1+4. J 


V1 , j+4  4  jj-ir—  (hi,j+l  ■  hi.j)  *  Vi,j+4 


Ayi »J+4 


(A6) 

(A7) 
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where  the  left  hand  side  of  (A6)  through  (A8)  are  the  unknowns  at  time  level  n+1 
and  u*,  V*,  and  H*  represent  all  known  quantities  at  time  levels  n  and  n-1,  cf. 

(Al)  through  (A3).  The  grid  index  in  the  x-direction  is  i,  in  the  y-direction,  j, 
with  integral  values  representing  the  h-grid.  Substituting  (A6)  and  (A7)  into  (A8) 
we  obtain 


11 /D  - 


i _  (h1»l.j'hi.j  .  h1J~hi-i.j\  a _ /N.j-H’N  ,j  _  ht ,j~hi  ,j-A 

(A9) 

H*‘  “  <uwu  -  *T,j.  J 


-  H?.o  - 


AX 


1  ,J 


1  ,J 


where  a=gHAt2.  At  the  western  boundary  i=h>,  and  Ui  is  known  at  the  necessary 
time  levels:  n-1,  n,  n+1.  Thus,  one  half  grid  distance  inside  this  boundary, 
instead  of  (A9) ,  we  obtain 


h  ,  .  i_  (  h2,J~hl. 
l.J  AX,  A  AX  3  A 

l»J 


h-  ‘',l'  1  ,  \i)  .  i 

gHAt/  Ay 

2*  ' 


-  HT,j  -  4t 


__  /.hu±r.,h.u 

I.M.J+h 


i . 


hl^:hLjrl) 


(A10) 


The  h-field  is  determined  in  a  similar  fashion  at  the  other  boundaries  and  the 
corners  of  the  h-grid.  When  the  boundary  conditions  are  included  in  the  forcing 
function  on  the  right  hand  side  of  the  equation,  the  problem  is  reduced  to  the 
solution  of  a  Helmholtz  equation  with  homogeneous  Neumann  boundary  conditions 
one-half  grid  distance  outside  the  boundary  of  the  h-grid.  Hence,  the  need  to  know 
the  true  value  of  h  or  its  normal  derivative  at  the  boundary  has  been  eliminated. 


73 


This  is  true  even  at  the  eastern  port  where  the  normal  flow  is  self-determined 


with  the  integral  constraint  that  the  total  transport  in  each  layer  exactly 
compensate  the  prescribed  inflow  through  the  southern  port  (see  Pig.  2) .  The 
normal  flow  at  the  boundary  of  the  domain  in  the  eastern  port  is  predicted  using 
the  x-component  of  the  momentum  equation,  (1) .  At  times  the  model  predicts  inflow 
for  part  of  the  lower  layer  at  the  eastern  port.  The  boundary  condition  used  for 
this  inflow  is  Ux=0,  eliminating  the  advective  term,  (Uu)Xr  and  setting 

Uxx  •  2  -  Vi)'"2  •  <A11> 


where  N  is  the  index  of  the  boundary  in  the  x-direction.  At  outflow  points  the 
computational  boundary  condition  used  for  the  normal  flow  is  Uxx*0*  In  addition  to 
eliminating  a  component  of  the  viscous  term  it  causes  Lilly's  (1965)  scheme  to  be 
exactly  equivalent  to  using  upstream  differencing  for  the  (Uu) x  term  at  outflow, 
when  we  perform  the  computation  at  time  level  at  t-At: 


(Uu)x  =  fUN(uN“uN_1)  +  uN^UN“UN-1^| /AX 


(A12) 


The  evolution  of  our  solutions  is  sufficiently  slow  that  gravity  waves  with 
significant  amplitude  are  not  excited,  so  the  reflection  of  these  waves  is  not  a 
problem.  The  tangential  flow  at  the  ports  is  set  to  zero  one-half  grid  distance 
outside  the  physical  domain  of  the  model.  This  weak  over specification  at  outflow 
eliminates  the  possibility  that  the  Loop  Current  will  exit  the  basin  at  unrealistic 
angles.  More  realistically,  outflow  should  occur  through  a  channel  modeling  the 
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Florida  Straits.  This  was  done  in  a  few  of  the  numerical  experiments. 

With  the  modifications  described  above,  we  have  sufficient  information  to 
calculate  all  the  terms  in  the  x-momentum  equation  required  to  predict  at 

the  eastern  port,  except  for  the  normal  pressure  gradient.  This  is  determined  by 
an  integral  constraint  on  the  total  transport  through  the  eastern  port  in  each 
layer.  When  the  integral  constraint  is  applied  in  the  determination  of  j  ,  the 
normal  pressure  gradient  is  assumed  to  be  uniform  across  the  port.  For 
approximately  geos trophic  flow  and  a  north-south  oriented  boundary  as  we  have  here 
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x  y 


fuxdy  +9hx  yr 


(A13) 


Hence,  a  uniform  presssure  gradient  is  a  useful  approximation,  if  the  integral  in 
(A13)  is  much  less  than  the  mean  pressure  gradient,  or  if  the  normal  pressure 
gradient  is  a  locally  anall  term  in  the  x-component  of  the  momentum  equation  at  the 
px>rt.  See  Hurlburt  and  Thompson  (1973,  1976)  and  Hurlburt  (1974)  for  discussion  of 
the  normal  pressure  gradient  in  an  east-west  oriented  cross-section. 

We  wish  to  calculate 
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where  u*  represents  all  the  other  quantities  in  the  x-momentum  equation  known  at 
time  levels  n,n-l.  Since  we  assure  hx^hx(y)  in  the  port 
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where  the  overbar  denotes  an  average  across  the  port.  In  our  case  we  must  express 
this  in  terms  of  transports.  Note  U=uh  +  u  V,  and  U*»u*h+u* 'h 1 ,  where  primed 
quantities  are  deviations  from  the  mean.  Since  uj|J+j  =  u^njn"^  +  constant,  u'«u*’, 

u  -  ui  =(U  -  U*)/h  (A16) 

and 


n+1  n-1  /n+1 

UN,j  =  hN-^,  j  vN 


U 


*n,n-l 

N.j 


(A17) 


an  equation  which  when  integrated  across  the  port  gives  the  desired  total  transport 

based  on  a  uniform  normal  pressure  gradient  across  the  port.  The  h  values  used  are 

one-half  grid  distance  into  the  interior  at  a  known  time  level,  and  ujj+1  is  the 

prescribed  total  transport  in  the  layer  through  the  port,  divided  by  the  width  of 

the  port.  This  transport  is  the  only  information  at  the  ports  required  at  a  future 

time  level  a  priori.  Since  it  is  prescribed  as  a  function  of  time,  we  can  predict 
n+1  n+1 

j  in  the  port  before  calculating  h^f j  using  the  Helmholtz  solver.  Thus,  the  solver 
is  provided  with  the  boundary  information  required  even  at  the  ports  where  the 
boundary  values  are  self-determined.  Once  we  have  h?+^  ,  the  semi- implicit 
calculation  of  U<+?  .  and  V<+1  ,  In  the  Interior  can  be  completed. 


The  behavior  of  the  open  boundary  condition  has  been  surprisingly  good.  When  a 
similar  boundary  condition  was  used  on  the  west  side  of  the  basin,  a  large 
amplitude  internal  Rossby  wave  in  a  westward  mean  flow  was  observed  to  pass  cleanly 
through  (unpublished  results) .  Although  ordinary  gravity  waves  are  reflected, 


Kelvin  waves  are  observed  to  pass  through  (Hurlburt,  Kindle,  and  O'Brien,  1976). 
Strong  jets  pass  through  (e.g.,  see  Section  4a,  Fig.  19a),  as  do  eddies  (unpublished 

results) .  In  no  case  in  the  several  applications  so  far  have  boundary  layers  with 


significant  amplitude  appeared  along  an  open  boundary  of  this  type.  Also  no 
smoothing  or  extra  damping  at  the  boundary  was  used  in  our  application. 

In  the  open  boundary  condition  with  self-determined  flow  the  mean  is 
constrained.  In  the  case  of  an  east-west  boundary,  deviations  from  the  mean  flow 
are  somewhat  constrained  by  the  tendency  to  form  a  Sverdrup  interior  (Hurlburt  and 
Thompson,  1973).  However,  this  is  not  true  for  a  north-south  oriented  boundary  as 
here,  allowing  the  possibility  of  unpnysical  growth  of  these  modes.  When  the 
entire  eastern  boundary  of  the  Gulf  of  Mexico  model  was  opened,  a  gyre  with  the 
north-south  scale  of  the  basin  developed  in  the  eastern  Gulf  on  a  time  scale  much 
longer  them  an  eddy  shedding  cycle  (see  Section  4a).  It  is  driven  by  the  strong 
Loop  Current  flowing  out  of  the  basin  along  the  southern  boundary,  a  physically 
reasonable  result.  When  the  model  domain  with  the  open  eastern  boundary  was 
extended  further  eastward,  the  growth  rate  of  the  gyre  was  increased  (Table  2, 

Cases  HG44,45),  perhaps  because  a  greater  percentage  of  the  perimeter  of  the  gyre 
within  the  model  domain  was  in  contact  with  the  strong  southern  boundary  current. 
This  demonstrates  an  unphysical  aspect  of  this  growth,  i.e.,  the  growth  rate 
dependence  on  the  percentage  of  the  gyre  within  the  model  domain. 

The  Helmholtz  equations  are  solved  using  a  vectorized  variant  of  Hockney's 
(1965)  method  by  Dan  Moore  (personal  communication)  for  homogeneous  Neumann 
boundary  conditions  one-half  grid  distance  outside  the  h-grid.  The  nonrectangular 
geometry  which  incorporated  the  Yucatan  and  Florida  Straits  was  included  using  a 
version  of  the  capacitance  matrix  technique  by  Dan  Moore  (personal  communication, 
see  Hockney  1968,  1970) .  This  requires  a  second  call  to  the  Helmholtz  solver.  The 
Yucatan  Peninsula  was  included  in  one  case,  in  which  the  Helmholtz  equation  was 
solved  using  a  version  of  the  stabilized  error  vector  propagation  (SEVP)  method  by 
David  Dietrich  (personal  communication) .  Another  version  of  this  method  is 
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described  by  Madala  (1978) .  Hie  implicit  part  of  the  finite  difference  equations 
for  the  gravity  waves  treat  the  gravity  wave  speed  as  a  constant  for  each  vertical 
mode.  Deviations  from  this  constant  occur  due  to  the  dynamics  and  the  topography. 
These  components  of  the  gravity  waves  are  treated  explicitly.  This  can  be  done 
without  introducing  conditional  stability,  if  the  deviations  always  act  to  slow  the 
gravity  waves  from  the  reference  value. 

This  layered  model  with  transports,  U,  V,  as  dependent  variables  handles 
topography  especially  well.  The  continuity  equation  is  linear  in  transport.  The 
topography  appears  multiplicatively  in  the  pressure  gradient.  The  topography  is 
differentiated  only  to  the  extent  that  it  affects  the  velocity  field  in  the 
advective  terms.  When  large  amplitude  topography  is  introduced,  restrictions  on 
the  time  step  and  the  eddy  coefficient  are  affected  only  to  the  extent  that  the 
topography  affects  the  amplitude  of  the  velocity  field  (advective  CFL  condition) 
and  the  grid  interval  Reynolds  number. 
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FIGURE  LEGENDS 


Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


Topography  of  the  22°C  isothermal  surface,  August  4  to  18,  1966 
(Alaminos  cruise  66-A-ll)  from  Leipper  (1970).  Subsequent  data 
(Elliot,  1979)  indicates  that  an  anticyclonic  eddy  separated 
from  the  Loop  Current  within  several  months. 

Bathymetry  of  the  Gulf  of  Mexico  based  on  U.  S.  Coast  and  Geo¬ 
detic  Survey  Chart  1007  and  soundings  on  file  at  the  Department 
of  Oceanography,  Texas  A&M  University.  From  Nowlin  (1972).  The 
rectangle  shows  the  approximate  domain  of  the  numerical  model. 

The  locations  of  the  inflow  and  outflow  ports  are  also  indicated. 
Bathymetry  of  the  idealized  Gulf  of  Mexico  model.  The  deepest 
water  is  at  3000  m  and  the  shallowest  topography  is  400  m  deep. 
The  contour  interval  is  ^0  m. 

Sequence  of  snapshots  of  the  PHA  (pycnocline  height  anomaly)  at 
70  day  intervals  showing  the  life  cycle  of  an  eddy  starting  at 
day  2210.  The  contour  interval  is  20  m.  In  all  the  figures 
dashed  contours  are  negative.  PHA  is  positive  downward.  The 
case  shown  here  uses  the  parameters  of  Table  1  and  the  topo¬ 
graphy  of  Fig.  3,  except  that  the  upper  layer  inflow  transport 
is  25  Sv,  the  lower  5  Sv. 

PHA  vs.  time  190  km  west  of  the  center  of  the  inflow  port  for 
3  cases:  (a)  case  shown  in  Fig.  4  which  includes  topography, 

(b)  standard  two-layer  flat  bottom  case  using  parameters  from 
Table  1,  (c)  standard  reduced  gravity  case  using  appropriate 
parameters  from  Table  1.  A  regular  progression  of  eddies 
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Fig.  6 


Fig.  7 


Fig.  8 


Fig.  9 


through  the  north-south  cross-section  is  shown  in  each  case. 

The  contour  interval  is  30  m. 

FSA  (free  surface  anomaly,  positive  upward)  at  5  cm  contour 
intervals  for  a  two-layer  flat  bottom  case  using  parameters 
from  Table  1  except  that  8=0  and  the  inflow  is  angled  27°  west 
of  normal.  The  solution  evolves  to  the  steady  state  shown  here  at 
day  1080. 

Snapshots  of  PHA  at  60  day  intervals  beginning  at  day  1320  for 
the  standard  reduced  gravity  case  after  B  was  reduced  to  zero 
with  a  time  constant  of  10  days  starting  at  day  1300.  The 
westward  spreading  of  the  Loop  Current  has  been  halted  and  the 
eddy  shedding  process  appears  frozen.  The  contour  interval  is 
20  m. 

PHA  for  (a)  the  mean  over  3  eddy  cycles  for  the  standard  re¬ 
duced  gravity  case,  (b)  the  deviations  from  the  mean  at  day 
1800,  (c)  a  steady  state  reduced  gravity  solution  (day  1800) 
obtained  by  increasing  the  eddy  viscosity  of  the  standard  case 
to  A  =  3  x  10  cm  sec  ,  and  (d)  the  steady  state  linear  solu¬ 
tion  (day  1080)  corresponding  to  the  standard  case,  obtained 
using  an  inflow  transport  of  .1  Sv.  The  contour  interval  is 
20  m  for  (a),  (b),  and  (c).  In  (d)  both  inflow  transport  and 
contour  interval  were  reduced  by  a  factor  of  200. 

FSA  for  4  baro tropic  flat  bottom  experiments.  The  appropriate 
parameters  from  Table  1  were  used  except  in  (a)  depth  =  200  m, 
inflow  =  10  Sv,  (b)  depth  =  300  m,  inflow  =  20  Sv,  (c)  depth  = 

300  m,  inflow  =  25  Sv,  (d)  depth  =  300  m,  inflow  =  30  Sv.  The 
contour  interval  is  5  cm.  Snapshots  are  at  (a)  day  720  (b)  day 
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710,  (c)  day  720,  (d)  day  1080.  In  panel  (a)  the  solution  is 
steady. 

Fig.  10  Snapshots  of  PHA  at  20  day  intervals  starting  at  day  1300  for 

6  2 

the  standard  two-layer  flat  bottom  case  except  A  =  3  x  10  cm  sec 
The  model  domain  has  been  extended  160  km  eastward  to  include  a 
representation  of  the  Florida  Straits.  The  contour  interval  is 
20  m. 

Fig.  11  Kinetic  energy  in  the  rectangular  part  of  the  basin  as  a  func¬ 
tion  of  time  for  each  layer  (upper  curve  for  upper  layer)  for 
3  cases:  (a)  standard  two-layer  flat  bottom  case  except 
A  =  3  x  10  cm  sec  and  the  domain  of  Fig.  10  is  used,  (b) 
standard  two-layer  flat  bottom  case,  (c)  standard  reduced 
gravity  case  except  A  =  3  x  10  cm  sec  .  In  case  (c)  the  upper 
curve  is  potential  energy.  The  value  of  IC  is  (a)  .7,  (b)  .5, 
(c)  5.5. 

Fig.  12  Panels  (a)  and  (b)  are  snapshots  of  PHA  at  extremes  of  an  eddy 
cycle  (days  1050  and  1080)  when  the  standard  two-iayer  flat 
bottom  model  has  inflow  27°  west  of  normal  and  the  lower  layer  i 
flow  is  reduced  to  zero.  Panel  (c)  is  PHA  vs  time  for  this 
case  at  a  longitude  190  km  west  of  the  center  of  the  inflow  port 
The  contour  interval  is  20  m. 

Fig.  13  Kinetic  energy  in  the  basin  vs  time  for  (a)  the  case  shown  in 
Fig.  12.  The  upper  curve  is  for  the  upper  layer  and  the  lower 
curve  for  the  lower  layer;  (b)  is  for  the  corresponding  reduced 
gravity  case.  The  value  of  IC  is  (a)  .55,  (b)  .9. 

Fig.  14  Potential  energy  spectral  densities  vs.  frequency  at  a  given 
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location  in  the  basin  for  4  cases:  (a)  the  standard  two-layer 

flat  bottom  case  at  x  =  1010  km,  y  =  366  km,  (b)  the  reduced 

6  2  -1 

gravity  case  with  A  =  3x10  cm  sec  at  x  =  1010  km,  y  =  366 

6  2  -1 

km,  (c)  the  two-layer  flat  bottom  case  with  A  =  3x10  cm  sec" 
at  x  =  1410  km,  y  =  291  km,  and  (d)  the  two-layer  flat  bottom  case 
with  inflow  27°  west  of  normal  confined  to  the  upper  layer  at 
x  =  1010  km,  y  =  234  km.  Spectra  were  computed  over  an  integral 
number  of  eddy  cycles  after  the  solutions  reached  statistical 
equilibrium. 

Fig.  15  North-south  section  of  PHA  in  meters  versus  time  190  km  west 
of  the  center  of  the  inflow  port  for  the  reduced  gravity  model 
with  inflow  centered  at  x  =1000  km  and  outflow  centered  at 
y  =  225  km.  The  contour  interval  is  30  m. 

Fig.  16  Snapshots  of  PHA  at  day  1080  for  (a)  the  standard  reduced  gravity 

case  and  (b)  the  case  shown  in  Fig.  15.  The  contour  interval  is  20  m. 

Fig.  17  Steady  state  PHA  (day  1080)  for  the  standard  reduced  gravity  model 
but  with  a  section  of  land  area  added  east  of  x  =  1440  km  and  north 
of  y  =  187.5  km.  The  contour  interval  is  20  m. 

Fig.  18  Regime  diagram  for  three  stability  regimes:  (E)  Eddy-shedding,  (W) 
steady  westward-spreading  (Fig.  8c,  d)  and  (N)  steady  source-sink 
(Figs.  6,  17)  as  a  function  of  Reynolds  number  (Re)  and  beta  Rossby 
number  (Rg).  Lower  case  letters  refer  to  barotropic  experiments; 
upper  case  letters  refer  to  reduced-gravity  experiments. 

Fig.  19  Snapshot  of  PHA  at  day  1080  (a)  and  y-t  phase  plot  of  PHA  190  km 
west  of  the  center  of  the  inflow  port  (b)  for  the  standard  reduced 
gravity  model  except  that  the  eastern  boundary  is  entirely  open  and 
the  east-west  extent  is  1760  km.  The  contour  interval  is  20  m  for 
(a)  and  30  m  for  (b). 
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Fig.  20 

Fig.  21 

Fig.  22 

Fig.  23 

Fig.  24 

Fig.  25 

Fig.  26 


Snapshots  of  PHA  at  day  1080  for  the  standard  reduced  gravity 
model  with  the  Yucatan  and  Florida  Straits  added;  (a)  for  stan¬ 
dard  outflow  port  width,  (b)  for  -halved  outflow  port  width  and 
(c)  for  doubled  outflow  port  width.  The  contour  interval  is  20 
m.  The  total  domain  size  is  1760  x  960  km. 

Energy  balances  for  the  three  models  of  Fig.  20  averaged  over 
an  eddy-cycle  and  over  the  rectangular  part  of  the  domain. 

Snapshot  of  PHA  at  day  1080  (a)  and  y-t  phase  plot  of  PHA  190 
km  west  of  the  center  of  the  inflow  port  (b)  for  the  standard  re¬ 
duced  gravity  model  with  the  Yucatan  Peninsula  represented  by 
a  rectangular  coastline  irregularity.  The  contour  interval 
is  20  m  for  (a)  and  30  m  for  (b). 

Steady-state  FSA  at  day  1080  for  (a)  the  flat  bottom  barotropic 
standard  experiment  and  (b)  the  same  experiment  with  bottom 
topography  of  Fig.  3  added.  Inflow  transport  was  30  Sv.  The 
contour  interval  is  2  cm. 

Snapshots  of  PHA  at  day  1980  for  (a)  the  standard  reduced 
gravity  model  but  with  25  Sv  inflow,  (b)  the  flat  bottom,  two- 
layer  model  with  parameters  identical  to  (a)  and  no  deep  water 
inflow,  and  (c)  identical  to  (b)  but  with  bottom  topography. 

The  contour  interval  is  20  m. 

Phase  plots  of  PHA  390  km  west  of  the  center  of  the  inflow  port 
for  the  three  models  shown  in  Fig.  24.  The  contour  interval  is 
30  m. 

Results  for  the  full  bottom  topography  experiment  with  inflow 
transport  T^  =  25  Sv  and  T^  =  0  until  year  6  when  T^  increases 
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to  10  Sv.  (a)  Domain-averaged  upper-layer  (top  curve)  and  lower- 
layer  kinetic  energy.  IC  =  1.5.  (b)  Steady  state  PHA  at  day 

2880.  The  contour  interval  is  20  m.  (c)  Lower-layer  pressure 

normalized  by  density  at  day  2880.  The  contour  interval  is 

2  -2 
.25  m  sec  . 

Fig.  27  (a)  Bottom  topography  of  the  idealized  Gulf  of  Mexico  when  the 

Campeche  Bank  and  the  bottom  slope  on  the  inflow  are  removed. 

The  deepest  water  is  3000  m  and  the  shallowest  water  is  400  m. 

The  contour  interval  is  ^00  m.  (b)  Domain-averaged  upper-layer 
(top  curve)  and  lower-layer  kinetic  energy  for  this  experiment. 

T.|  =  25  Sv  and  T2  was  increased  from  0  to  10  Sv  after  3  years. 

IC  =  .95. 

Fig.  28  Steady-state  results  (day  1080)  for  the  experiment  shown  in  Fig. 

27.  (a)  Bottom  topography  (dashed)  contoured  at  250  m  intervals 

and  depth  of  the  interface  (solid)  contoured  at  20  m  intervals 
in  a  500  x  500  km  section  in  the  southeast  corner  of  the  model, 
(b)  Contours  of  V  .  v-j  •  (c)  Contours  of  v^  .  V  h^ .  (d) 

Contours  of  (g/f)  J  ^.n^).  The  contour  interval  in  b  -  d  is 

-4  -1 

30  (x  10  cm  sec  ). 

Fig.  29  Steady-state  (day  1080)  mass  transport  vorticity  term  balances 
for  the  experiment  shown  in  Fig.  27  at  y  =  168.75  km  across 
the  Loop  Current  from  x  =  1000  km  to  x  =  1600  km  for  (a) 
upper  layer,  (b)  lower-layer,  and  (c)  the  sum  of  (a)  and  (b). 

From  equations  (21-23):  diffusion  (solid),  planetary  vorticity 
advection  (long  dashes),  pressure  torques  (short  dashes),  and 
nonlinear  terms  (dots  on  solid). 
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Fig.  30  Potential  energy  spectral  densitites  versus  frequency  for  the 


standard  reduced  gravity  model  at  x  =  810  km,  y  =  421.9  km 
with  (a)  steady  forcing  and,  superimposed  on  the  mean  (b) 
an  annual  forcing,  (c)  an  eight-month  period  forcing,  and  (d) 
a  semi-annual  period  forcing.  The  mean  transport  was  20  Sv; 
the  sinusoidal ly-varying  amplitude  was  10  Sv  peak  to  peak. 
Spectra  were  computed  over  an  integral  number  of  eddy  cycles 
after  the  solutions  reached  statistical  equilibrium. 

Fig.  31  Phase  plots  of  PHA  390  km  west  of  the  center  of  the  inflow 

port  for  (a)  the  eight-month  period  forcing  experiment  and  (b) 
the  semi-annual  period  forcing  experiment.  The  contour  interval 
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(b)  Reduced  Gravity  -  Straits  -  Half  Width  Outflow 
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Current.  Typically  the  models  were  integrated  3  to  5  years  to  statistical 
equilibrium  on  a  1600  x  900  km  rectangular  domain  with  a  resolution  of 
20  x  18.75  km.  Prescribed  inflow  through  the  model  Yucatan  Channel  was  com¬ 
pensated  by  outflow  through  the  Florida  Straits. 

A  long-standing  hypothesis  is  that  the  Loop  Current  sheds  eddies  in  re¬ 
sponse  to  quasi -annual  variations  in  the  inflow.  We  find  that  the  Loop  Current 
can  penetrate  into  the  Gulf,  bend  westward,  and  shed  realistic  anticyclonic 
eddies  at  almost  an  annual  frequency  with  no  time  variation  in  the  inflow.  In 
this  regime, the  eddy  shedding  rate  depends  on  the  internal  Rossby  wave  speed, 
an  eddy  diameter  derived  from  conservation  of  potential  vorticity  on  a  6-plane, 
the  angle  of  the  inflow,  and  to  a  lesser  extent  on  the  Reynolds  number.  Eddy 
shedding  can  be  prevented  by  reducing  the  Reynolds  number  sufficiently.  How¬ 
ever,  the  Loop  Current  still  spreads  far  westward.  The  steady  state  solution 
for  a  highly  viscous  case  was  found  to  be  almost  the  same  as  the  mean  over  an 
eddy  cycle  for  a  lower  viscosity  case  which  shed  discrete  eddies  of  large  am¬ 
plitude.  Eddy  shedding  and  westward  spreading  of  the  Loop  can  be  prevented 
at  higher  Reynolds  numbers  when  the  beta  Rossby  number,  Rg=vc/(£LpO  -2,  where 
the  appropriate  length  scale,  Lp,  is  one-half  the  port  separation  distance  and 
v^  is  the  velocity  at  the  core  of  the  current.  Differential  rotation  (f)  is 
also  of  great  importance  in  determining  the  diameter  and  westward  speed  of  the 
eddies  and  the  penetration  of  the  Loop  Current  into  the  Gulf.  In  a  few  of  the 
two-layer  experiments ,  baroclinic  and  mixed  instabilities  were  encountered, 
but  experiments  dominated  by  a  horizontal  shear  instability  of  the  internal 
mode  produced  the  most  realistic  results.  For  sufficiently  high  Reynolds 
numbers  the  shear  instability  occurred  in  both  the  barotropic  and  reduced 
gravity  models.  However,  for  realistic  parameter  values  eddy  shedding  occurred 
in  the  two-layer  and  reduced  gravity  models,  but  not  in  the  barotropic  model. 

Consistent  with  potential  vorticity  conservation  dynamics,  the  Loop 
Current  and  its  eddy  shedding  behavior  was  quite  insensitive  to  the  location 
and  width  of  the  inflow  and  outflow  ports,  so  long  as  the  western  boundary  did 
not  interfere  with  the  shedding  process  and  the  ports  were  not  separated  by 
much  less  than  W, '2  times  the  theoretical  eddy  diameter,  i.e.,  when  Rg<2.  Whenj 
the  entire  eastern  boundary  was  opened,  the  outflow  remained  confined  to  a  cur¬ 
rent  adjacent  to  the  southern  boundary.  Also,  while  the  solution  depends  on 
the  maximum  velocity  at  inflow,  it  is  relatively  insensitive  to  the  shape  of 
tne  inflow  profile. 

In  the  presence  of  significant  deep-water  inflow  through  the  Yucatan 
Straits,  bottom  topography  may  prevent  Loop  Current  penetration,  westward 
spreading,  and  eddy  shedding.  In  these  cases  the  interaction  between  the 
bottom  topography  and  the  pressure  field  near  the  Florida  Shelf  results  in  a 
near  balance  between  the  pressure  torques  and  the  nonlinear  terms  in  the  mass 
transport  vorticity  equation.  When  the  Yucatan  Straits  deep  water  inflow  is 
reduced  or  the  Florida  Shelf  is  moved  to  the  east,  the  eddy  shedding  reappears. 
A  kinematic  analysis  shows  that  a  sufficiently  strong  current  following  f/h 
contours  of  the  Florida  Shelf  and  intersecting  the  Loop  Current  at  large 
angles  can  locally  prevent  northward  penetration  of  the  Loop  Current  and 
effectively  reduce  the  port  separation.  Thus,  the  effect  of  the  Florida  Shelf 
is  similar  to  cases  in  the  reduced  gravity  model  where  the  ports  are  too  close 
for  eddy  shedding  to  occur,  i.e.,  when  Rg  2.  Bottom  topography  also  inhibits 
development  of  baroclinic  instability,  yielding  solutions  more  closely  re¬ 
sembling  those  from  the  reduced  gravity  model  than  from  the  two-layer  flat 
bottom  model.  However,  movement  of  the  shed  eddies  is  significantly  modified 
by  the  introduction  of  topography. 

In  the  presence  of  realistic  time  variations  in  the  upper  layer  inflow, 
the  eddy  shedding  period  is  dominated  by  the  natural  period  rather  than  the 
forcing  period,  although  the  influence  of  the  latter  is  not  negligible. 
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